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This study explores the optimization method of simulated annealing for use in an 
analytic tool for counter-narcotics analysis. A model is developed, employing RAND 
Corporation’s Simulation of Adaptive Response model as an objective function evaluator, 
which optimizes interdiction asset locations relative to a Sample smuggling network in the 
Caribbean region. In addition to asset location optimization, the response of the model 
to changing numbers of assets is also tested. Results indicate that this methodology has 
potential for use in the counter-narcolics program, and perhaps other network interdiction 
applications. Further research and testing are recommended. Military drug interdiction 
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I. INTRODUCTION 


A. BACKGROUND 

The Defense Authorization Act of 1989 designated the Department of Defense 
(DoD) “as the single lead agency of the Federal Government for the detection and 
monitoring of aerial and maritime transit of illegal drugs into the United States." The Act 
also tasks DoD with expanded support to civilian Law Enforcement Agencies (LEA’s) 
involved in drug interdiction. According to Secretary of Defense Dick Cheney, 

The Department of Defense is an enthusiastic participant in the nation’s drug 
control effort. We have significant resources at our disposal. We can make a 
substantial contribution to our national effort if we use our assets intelligently and 
efficiently. [Ref. 1: p. 4] 
To these ends, the President’s National Drug Control Strategy for 1991 estimated the 
DoD budget for drug interdiction at nearly $1.1 billion, more than double the actual 
expenditures for 1989. Obviously, the potential contributions by DoD are being taken 
very Seriously. 

In 1988, the RAND Corporation developed the Simulation of Adaptive Response 
(SOAR) model to analyze the effect of military involvement in drug interdiction. SOAR 
is a Stochastic network simulation of drug smuggling involving theoretical routes for drug 
shipments and the response of smugglers to being interdicted. Expanded involvement of 


the military in drug interdiction was evaluated by examining the effect that increased 


interdiction probability on the theoretical routes might have on the import price of illegal 


drugs. The conclusions of the authors of the RAND study (Gordon Crawford, Peter 
Reuter, et al.) were not heartening. Although the limitations of the SOAR model caused 
it to be biased toward successful interdiction, it was found that “interdiction must be very 
stringent indeed to. greatly affect U.S. drug consumption." [Ref. 2: p. 73] 
However, the model is simplistic in its approach to smuggling routes and the presence of 
interdiction forces. Routes are not associated with actual geographic areas due to a lack 
of available data, while interdiction assets are represented solely by a probability of 
interdiction on the generic routes. This simplicity leaves room for expansion of the model 


and further analysis. 


B. PROBLEM ANALYSIS 

Government-sponsored studies prior to that undertaken by RAND did not 
concentrate on the adaptation of smugglers, nor explicitly examine the potential for 
military contribution to the drug interdicuon program [Rel. 3: p.4]. As part 
of a comprehensive study, SOAR was designed to aid in a policy analysis of the overall 
interdiction program. 

As a policy tool, SOAR did not require a geographic structure or specific 
representation of interdiction assets. For this reason, SOAR is not useful in operational 
planning or analysis. RAND’s conclusions, at a policy level, indicate that military 
involvement in drug interdiction may have limited effectiveness [Ref. 3: p. 130). 
However, the Department of Defense has a major role in the program, and would benefit 


from the availability of an analytic tool which could be used to evaluate the effects of 


operational planning and to explore, with a consistent measure, the effectiveness of 
different force and capability alternatives. 

Interdiction problems are frequently modeled as flow networks for which 
Optimization methods are used to determine which arcs or nodes should be cut or 
removed to provide an optimum result. In the case of a smuggling network, transit routes 
may be represented by arcs which are not well defined, or are completely unknown to the 
interdictors, making it difficult to entirely sever the arc. Even if it were possible, the 
smugglers themselves will eventually realize that the arc has been cut and will move to 
another route, adapting their methods to the interdiction. The SOAR model simulates this 
type of response on the part of drug smugglers, making it a fairly obvious choice to 
measure the effectiveness of interdiction on a smuggling network. 

Because routes used by smugglers may not be clearly defined, it is unrealistic to 
merely place assets on arcs representing segments of the routes. By developing a 
geographic structure for the routes, assets can be located by latitude and longitude, and 
the effects of assets on the routes (1.e., the probability of interdiction on routes) can be 
measured as a function of their range from the routes. This also would allow assets to 
affect multiple routes, for example, routes passing through a geographically restricted 
area, or chokepoint. Finding the optimum locations for assets relative to a smuggling 
network with stochastic flows will be difficult. As the locations of assets change, flows 
through the routes will also change, resulting in local and global optima which vary with 
the locations of assets. There are also likely to be many such optima, a result of different 


assets contributing the same probability of interdiction to routes or being located in 


positions previously tested by other assets. In addition, the number of possible locations 
for assets is infinite, requiring some control over the locations tested if discrete 
Optimization is to be used. This control must include geographically feasible areas, but 
should also reflect operational constraints, such as a requirement for certain assets to 
Operate within a set range from each other or a base, and political constraints, such as 
foreign airspace or waters restnctions. 

An optimization model with the properties described in this analysis, which retains 
enough realism to be considered valid, would be of great benefit at the operational level 


of the military drug interdiction program. 


C. METHODOLOGY 

We will propose a solution to the need for an analytical tool which may be used for 
analysis and planning at the operational level of the military drug interdiction program. 
An exploration of the use of simulated annealing to optimize the location of counter- 
narcotics assets with respect to a drug smuggling network is presented. The objective 
function is an outcome of the modified SOAR model, in which the routes used by the 
smugglers are modeled with a geographical structure, and the physical location of 
interdiction assets (eg. ships, aircraft, aerostat balloons, observation posts, etc.) are used 
to determine the probability of interdiction in an area. From the SOAR model we extract 
the total cost to smugglers, a measure of effectiveness for one major goal of the overall 
drug interdiction program, reduction in consumption of illegal drugs. Interdiction raises 


the risks involved with drug smuggling, thereby increasing the costs associated with 


personnel pay, equipment and drug replacement, etc. Eventually, this results in higher 
consumer prices for drugs. Economically, higher prices should cause a decrease in 
demand for drugs in the long run, that is, lower consumption by Americans. Increased 
crime among addicts to compensate for higher prices is not addressed. 

Use of SOAR as the objective function evaluator in a model which exhibits 
numerous local and global optimal asset locations eliminates the potential use of gradient- 
Optimization methods. It is for this reason that we chose to apply simulated annealing 


as the optmization method to maximize the total cost to the smugglers. 


D. SCOPE 

This thesis briefly discusses the flow of cocaine from South and Central America 
into the southern United States, and the role the Department of Defense plays in 
interdiction efforts within this geographical area. 

The model presented is capable of evaluating the use of assets in the Atlantic, 
Pacific and Southwest Border regions; however, for the numerical results presented in this 
thesis, only routes and assets in the Caribbean and Gulf of Mexico are examined. Much 
of the information that is available about known smuggling routes and the assets involved 
in counter-narcotics operations is highly classified. In the interest of maintaining an 
unclassified example, some data used in the model is taken from the original SOAR 
model, while information on smuggling routes and counter-narcotics assets is intended to 
be representative only, and should not be construed as data representing actual 


circumstances. 


Il. THREAT AND INTERDICTION EFFORTS OVERVIEW 


A. CURRENT THREAT OF ILLEGAL DRUG TRAFFICKING 

A number of indicators show that progress has been made in the war against drugs. 
Overall drug use, drug-related medical emergencies, and student attitudes toward drug use 
have all declined past the goals set in the 1989 National Drug Control Strategy, shown 


in Table 1. [Ref. 4: p. 5] 


Table 1 NATIONAL DRUG CONTROL STRATEGY NATIONAL OBJECTIVES: Statistical 


significance ranges from | to 30%. 
_ com. | _ ACTUAL 






FREQUENT COCAINE USE 50% - 23% 
reduction in actual decline 
rate of increase 






























CURRENT ADOLESCENT COCAINE USE - 20% - 49% 
DRUG-RELATED MEDICAL EMERGENCIES 
STUDENT PaTuDEs Towa 





Cocaine-specific indicators showed decreased purity and significantly higher prices 
in many areas of the country for the past year. Also, since 1985, "... reported current 
cocaine use (use in the past month) decreased some 72 _ percent." 


[RetoS: peat 


The effect on smugglers themselves may be deduced from the surrender of a 
number of Colombian drug lords and the level of effort smugglers put forth to get their 
product into the United States. Examples of this effort include the tunnel between 
Mexico and Arizona discovered in May 1990 [Ref. 6: p. Al] and the 
occasionally high tech hiding places used in vehicles, boats, and aircraft. 

On the other hand, many problems still exist which illuminate the threat to national 
security declared by President Reagan in 1986. Drug-related violence and murders remain 
at high levels within the United States, as does the amount of cocaine being produced 
throughout South America. Table 2 displays production levels for the three major 
producing countries from 1988 to 1990, the last year for which data was available. 


Table 2 COCAINE PRODUCTION ESTIMATE: Potential Cocaine HCI Production by Country 
(tons). Sources: NNICC Reports, 1989 and 1990. 


78 1-996 771-980 





Smugglers’ high tech methods display their increasing adaptability to interdiction 
efforts, as well as the effects of interdiction. Some smugglers are showing that they are 
willing (or crazy enough) to take on interdiction forces, as evidenced by incidents along 


the border with Mexico where smugglers have exchanged small arms fire with U.S. 


military forces [Ref. 7: p. 1], and the attempt in April 1990 by smugglers 
linked to the Medellin cartel to buy Stinger missiles [Ref. 8: p. 3]. 

The evolution of narco-terrorism may be proof enough that we are indeed facing 
a war which must be fought by our armed forces outside the United States. Many 
revolutionary/terrorist organizations, such as M-19 and FARC in Colombia and Shining 
Path in Peru, are intimately linked to the drug trade, forming the basis for this narco- 
terrorism. These groups have shown decided anti-American leanings, and have used 
violence and threats against high level U.S. officials. One-time M-19 leader Ivan Mariano 
Ospina said in December 1984, "May these threats be carried out and may they be carried 
out in the entire world against the rapacious imperialism that lives at the cost of misery 
of _ exploited. = pecple]. SalReie2]. Chao) Even state leaders have 
been involved, most recently illuminated by General Noriega in Panama, but also 
including Fidel Castro who has acknowledged the use of cocaine as a weapon, "We are 
going to make the people up there (the United States) white, white with cocaine [Ref. 9: 
oe SVAN 

Clearly, there are some signs of success in the drive to reduce illegal drug 
trafficking and abuse. However, it is also plain that much more effort is required in all 
phases of the war on drugs before the situation can be considered under control, let alone 
resolved. Involvement of the Department of Defense in reducing the flow of drugs into 
this country 1S only one part of the overall struggle, but, as Secretary Cheney said, the 


military can have a great impact in this effort. 


B. WHAT IS DOD’S INVOLVEMENT NOW? 

Since the Department of Defense’s first, somewhat reluctant, involvement in 
counter-narcotics in the 1970's, the military has been called upon for increasing support 
each year. The Defense Authorization Act of 1989 presented DoD with a primary role 
in the drug war. We now see all segments of the military accepting counter-narcotics 
assignments, including operational missions, such as ship and aircraft patrols with 
embarked law enforcement personnel, intelligence gathering and surveillance, and material 
support. The overall DoD chain of command for counter-narcotics is shown 


in Figure 1. 


SECDEF 
JCS: CNOD 


CINCLANT CINCPAC CINCSOUTH CINCFOR NORAD 
JTF4 JTFS5 JTF6 





Figure 1 DoD Counter-narcotics Chain of Command. 


Obviously there is high-level attention being paid to the problem. The three Joint 
Task Forces command the military assets involved in actual day-to-day interdiction 


operations. 


1. Joint Task Force 4 

Located in Key West, Florida, JTF4 is in charge of operations in the Atlantic 
area of responsibility (AOR). As the first Joint Task Force, their mission included the 
creation of a joint fusion center for tactical intelligence and the communications necessary 
to assemble the data and provide a finished product to all users. JTF4 conducts their own 
detection and monitoring operations in the Atlantic AOR, and coordinates the operations 
of other agencies, such as the Coast Guard, Customs Service and local law enforcement 
activities. 

Although JTF4 has no dedicated assets, they do have tactical control of 
various personnel, ships, patrol and intercept aircraft, and intelligence collecting assets 
assigned from other DoD commands, the Coast Guard and Customs Service. 


[Ref. 10: p. 78] 


2. Joint Task Force 5 
Located in Alameda, California, JTF5 has a similar responsibility and mission 
in the Pacific AOR as that of JTF4 in the Atlantic. They conduct organic operations to 
detect and monitor aircraft and surface vessels suspected of smuggling, as well as 
coordinating other agencies’ operations within their AOR. JTF5’s mission includes 


integrating into the intelligence communications network and providing drug-related 
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intelligence to law enforcement agencies as appropriate. [Ref. 11] As is the case 
with JTF4, JTF5 has no dedicated assets, but takes tactical control of various assets 


assigned from DoD commands, the Coast Guard and Customs Service. 


3. Joint Task Force 6 

Located at Ft. Bliss, Texas, JTF6 has a much different mission from that of 
JTF’s 4 and 5. Specifically, "Joint Task Force Six plans and coordinates all DoD support 
requested by Federal, State and local law enforcement agencies within the Southwest 
Border region.” [Ref. 12] 

Requests for support usually come through Operation Alliance, an agency 
comprised of representatives from Federal, State, and local law enforcement agencies of 
the Southwest Border states: Texas, New Mexico, Arizona, and California. Support 
provided includes observation posts, reconnaissance, terrain denial operations designed 
to deter smugglers from using areas along the border, training teams, military drug 
detection dogs, and construction. Personnel from all military services and various 


National Guard units have been called upon to perform these functions. 


C. WHERE MAY DOD’S INVOLVEMENT LEAD? 

There are as many different opinions concerning DoD’s future involvement in drug 
interdiction as there are people contemplating it. These opinions range from complete 
disassociation to complete involvement in all aspects, not just detection and monitoring. 
Many rational, and some not So rational, arguments are given for both extremes and 


everything in between. 


1] 


One of the main arguments against the use of military forces in counter-narcotics 
is the Posse Comitatus Act (Title 18 U.S. Code, Section 1385), which prohibits the direct 
use of military forces for civil law enforcement activities within the United States. 
Because of this law, military forces involved in drug interdiction operations may not 
perform search and seizure or make arrests. The military is limited to its monitoring and 
detection mission and support of law enforcement activities. 

Another strong argument is that drug interdiction operations have a detrimental 
effect on DoD’s primary mission to maintain a military capable of fighting a war to 
defend the U.S. or its interests. This argument contends that the fight against drugs is not 
truly a war, and that performing area patrols and ship boardings, detection and intercept 
of slow moving aircraft, and patrolling a border area are not suitable training operations 
for the military. Following the end of the Cold War (1.e., the reduction of the Soviet 
threat) drug smuggling takes its place among "principle threats ... that replace the old 
Soviet threat [Reta This means that drug interdiction 
Operations need no longer be construed merely as training evolutions, but are actually part 
of the military's primary mission. 

A third argument is the cost of military drug interdiction operations. It was 
estimated that the Air Force spent approximately $433,000 per drug bust in 1986. 
[Ref. 14: p. 7] Many believe that the $1 billion-plus military dmg 
interdiction budget could be better spent on education and treatment programs. These 
programs are extremely important, since interdiction alone cannot halt the problem of 


drug abuse. However, education and treatment are more effective when drugs are more 
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expensive and less readily available. "Most leaders of prevention and treatment programs 
recognize this; their task is made easier when drug enforcement works." 
[Ref. 15: p. 3] The goal of interdiction is to “complement and support our 
international drug control activities and domestic law enforcement programs, ... to create 
an integrated supply reduction program." [Ref. 15: p. 65] Interdiction functions not only 
to seize drugs destined for the U.S., but to act as a deterrent to potential drug smugglers, 
and to force the price of illegal drugs to nse through increased risks the smugglers face. 
In addition, "Interdiction provides a highly visible sign to other nations of U.S. interest 
in reducing drug use." [Ref. 3: p. 2] This is a typical extension of the age-old Navy 
mission of showing the flag. Thus, the cost per seizure of military interdiction operations 
does not encompass the whole military contribution to the overall drug control strategy. 

Plainly, the military has an important role to play in combatting the drug problem 
faced by the United States. While no one wants to see the military take on civilian law 
enforcement roles, DoD has the manpower, assets, technology, and money to assist and 
supplement Federal, State, and local law enforcement agencies in reducing the supply of 
drugs. 

While DoD maintains its role in detection and monitoring of illegal drugs entering 
the country, there are many who are pushing for an even greater military involvement. 
In 1990, Senator Mitch McConnel (R-KY) proposed a plan for shooting down suspected 
airborne smugglers who refuse to land for inspection [Ref. 16: p. Al9]. 
Although this plan was not put into effect, it did have a large amount of support. 


Military and LEA personnel are providing training and technical assistance for anti-drug 
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and counterterrorist forces in a number of production and transit countries. Operation 
Bahamas, the Turks and Caicos Islands (OPBAT), a multinational strike force in place 
since 1983, has been very successful in preventing smuggling flights into the Bahamas. 
It is certainly conceivable that similar strike forces could be established with other 
countries. 

The completion of the aerostat balloon network and the Canbbean Basin Radar 
Network, as well as proposed over-the-horizon radar systems, would significantly increase 
the ability to detect and monitor air traffic in smuggling areas, and act as a deterrent to 
smugglers. Aerostats are tethered radar-equipped balloons, flown at about 10,000 feet, 
which "provide surveillance of oceanic and land areas right down to the surface." 
[Ref. 17: p. 7], with a range of approximately 200 nautical miles. The 
Caribbean Basin Radar Network is a series of ground based radars designed to provide 
nearly complete coverage of the Caribbean basin. Figure 2 displays the proposed 
coverage of the aerostat network [Ref. 18: p. 84]. 

RADM Walter Leland, USCG, said, "We know that when the aerostats are flying 
the traffic stops." [Ref. 19: p. 7] However, there are widely differing 
reports on the operational reliability of the aerostat balloons. Some Air Force and 
Customs Service officials call the system too expensive and fragile, claiming the balloons 
are not operational more than 60% of the time because of Severe weather and 
maintenance problems _ [Ref. 20: p. 15]. At the same time, Pentagon 
figures showed the balloons to be operational 70-75% of the time. 


[Ref op sia Such discrepancies cast doubt on the effectiveness of 
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Figure 2 Aerostat Radar Network. 


these systems, in turn making it harder to get funding to continue and expand the 
program. 

Finally, there have been a number of popular novels recently published which 
provide more speculation about the possibilities of the military's future role in drug 
interdiction. Most notable of these are Tom Clancy's Clear and Present Danger and Dale 
Brown’s Hammerheads. 

Although many widely ranging proposals for military involvement in counter- 
narcotics and numerous arguments against any use of the armed forces in drug 
interdiction exist, with the current trend toward downsizing the armed forces it is difficult, 


if not impossible, to predict exactly what DoD’s future role in the nation’s drug control 


I5 


Strategy will be. It is fairly certain however, that the U.S. military will continue to be 
involved. A recent report by Rep. Les Aspin (D-Wis.), Chairman of the House Armed 
Services Committee, foresees the continuation of a military role: 


The drug problem must be attacked from both sides -- both the demand for drugs 
and the supply of drugs must be reduced.... Although the United States will 
continue to rely primarily on nonmilitary means in its effort to curtail the 
international trade in drugs, Amenicans will also want U.S. military forces engaged 
in the struggle with drug traffickers, primarily in an interdiction role. 

[Ref. 13: pou 
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lil. SIMULATION MODEL 


A. THE ORIGINAL SOAR MODEL. 


1. Background. 

The Simulation of Adaptive Response (SOAR) model is part of a RAND 
Corporation study on the drug interdiction program, undertaken at the request of the 
Under Secretary of Defense for Policy in February 1986. This study culminated with the 
publication of analysis results in three RAND reports. 

The primary conclusion of RAND’s effort was that: 

Increasing drug interdiction efforts are not likely to greatly affect the availability 
of cocaine in the United States....This conclusion is driven primarily by the small 
share of total drug distribution costs that are accounted for by the smuggling sector. 
Only about 10 percent of the final price of cocaine comes from smuggling costs and 
profits. [Ref. 3: p. x1] 
This conclusion does not take into account other reasons for military involvement in drug 
interdiction, such as those discussed in the previous chapter. It does, however, suggest 


that the need exists for a tool to be used in counter-narcotics analysis. RAND’s SOAR 
model provides a reasonable starting point for developing such a tool. 
2. Model Description. 
SOAR is a dynamic network simulation which allows drug smugglers a 
number of air, sea, and land smuggling routes. These routes are generic in that they are 


not associated with particular geographic routes. The route used for a particular shipment 
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of drugs is probabilistically selected, depending on the costs of shipping drugs via that 
route and the expected risk costs for shipping via the method associated with that route 
based on a perceived probability of the shipment being seized. This probability is 
estimated by a time-weighted history of unsuccessful shipments on that route, and will 
change with time as shipments succeed or are interdicted. 

Time between shipments is an exponential random variable with an input 
mean. The quantity of each shipment is an input constant, aS are the capacity and costs 
of the routes and shipping methods. SOAR does not attempt to model the adaptation of 
changing shipment sizes in response to interdiction severity. The capacity of the route 
and the size of the shipment are used to determine the number of trips that each shipment 
will require. These trips occur on the same route at the same time, leading to a Saturation 
factor which increases the probability that the shipment 1s seized. 

Interdiction in SOAR is modelled by a general probability of interdiction for 
each route. "Military resources are treated...simply as means for augmenting particular 
activities, thereby raising the probability of success in those activities." [Ref. 3: p. 8] 
This probability is an input value and can be changed for predetermined time periods 
(phases). When multiple trips are required, the probability of interdiction on the route, 
for that shipment, is increased as a result of the saturation factor. 

The smugglers’ overall goal 1s to ship as large a quantity of drugs as possible 
at their perceived lowest cost. Costs to the smuggler include the cost of the drug, 
operating costs for the method of shipment, the cost of replacing lost assets if a shipment 


is seized, and pay to personnel. Smuggling personnel are assumed to be risk averse, 
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therefore pay is assumed to change as the square of the perceived risk involved in using 
a particular route. Increasing interdiction rates will affect the overall cost to the 
smugglers by raising each of these costs. 

For their model, the RAND authors attempted to estimate: 

¢ Quantities shipped, by route, in a given year. 
¢ Number of shipments, by route, in a given year. 
¢ Number of vessels identified as suspicious, by route, in a given year. Of those, the 
number pursued, by route in a given year. Of those, the number resulting in 
seizures, by route, in a given year. 
¢ Estimates of the compensation resulting from the likelihood of prison. 
¢ Estimates of smugglers’ nonrisk compensation and profits. [Ref. 2: p. 9] 
As can be expected, little concrete data was available, resulting in the use of informed 
guesses for many of the input parameters. 

Output of the model includes an echo of the input parameters and the statistics 
described in Appendix A, Sections A.2. and B.1. respectively. The output is an average 
of an input number of runs desired. Ten runs were used in the RAND analysis, providing 
"an adequately precise estimate of the overall means, especially compared with the 
imprecision of some of the input data and some of the assumptions incorporated in the 
model." {Ref. 2: p. 24] This output was not used in our analysis, but is included in the 
appendix for information. 

Computational details of the SOAR model are described in Chapter IV of 


Simulation of Adaptive Response, and are reasonably easy to interpret from the model 


itself. 
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B. THESOAR SIMULATION AS THE OBJECTIVE FUNCTION EVALUATOR. 


1. Measures of Effectiveness. 

In evaluating drug interdiction efforts, many measures of effectiveness have 
been proposed. Civilian law enforcement agencies are concerned with numbers of 
seizures and smugglers prosecuted. News reports regularly include the street value of 
drugs seized. Because of its unique role of detection and monitoring and support to law 
enforcement, and the fact that the armed forces are prohibited from direct participation 
in the arrest and seizure phases of interdiction, the military's measure of effectiveness 
(MOE) is based on the number of supporting missions provided. This MOE is not very 
useful for a model which attempts to find optimal locations for interdiction assets. 
Although it indicates the level of military support to reducing the flow of drugs into the 
United States, it provides no quantitative measure of the effectiveness of the overall drug 
interdiction program. 

The quantity of drugs seized is a popular MOE, as it provides an evaluation 
of how much drugs are prevented from entering the country. However, it is possible that, 
when interdiction efforts are most successful, no drugs at all would be seized. That is, 
shipments could be deterred, smugglers could be forced to use the most expensive 
methods of shipping to avoid interdiction, etc. 

The major goal of drug interdiction is to reduce the consumption of illegal 
drugs. Military assets are an integral part of increasing the effectiveness of interdiction, 
thereby helping to achieve this goal. In addition to seizing drugs and equipment, 


interdiction increases the risks which smugglers face and causes them to change their 
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methods of operation. All of these effects of interdiction result in higher costs to the 
smugglers, costs which are eventually passed on to the consumers of illegal drugs. 
Although some users are likely to be insensitive to increased drug prices, "the body ot 
economic support for a relationship between price and consumption is too strong to allow 
much doubt that, at least in the long run, higher prices will lead to lower consumption." 
[Ref. 3: p. 20] Therefore, if it can be measured or estimated with some confidence, cost 
to the smugglers is an MOE which reflects the overall goal of the drug interdiction 
program. 

As computed by the SOAR model, total cost to the smugglers includes the 
costs of unsuccessful shipments, thereby incorporating the MOE of the amount of drugs 
seized. Parameters for an alternate objective function value could easily be extracted 
from the simulation, or the optimization algorithm could be applied to any desired 
objective function. 

The simulation used as the objective function evaluator for the opumization 
model presented in this thesis 1s a slightly modified version of the onginal SOAR model. 
The changes involve the attempt to model smuggling routes with a geographical structure 


and calculating the probability of interdiction based on the location of military assets. 


2. Model Assumptions. 
In any model, many simplifying assumptions are required. This is especially 
true in attempting to model a process of this scale. Perhaps the biggest assumption lies 
in the use of the SOAR model outcome as the objective function to be optimized. That 


is, we assume that the simulation provides a reasonable estimate of the total cost to the 
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smugglers. This includes accepting as reasonable the assumptions made in constructing 
SOAR, which are described in the RAND publications. 

Another major assumption involves modelling the drug smuggling routes. 
Except for the natural chokepoints in the Caribbean, little is known about the exact routes 
being used. Therefore, the routes presented in this model are assumed to be 
representative of area routes of unknown width. In the Mona Passage, Yucatan Channel 
and other chokepoints, the width of sea routes 1s constrained by the presence of land on 
either side. However, in the Pacific, Southwest Border area, and for many air routes, the 
width could be hundreds of miles. Assets which provide coverage of a route are assumed 
to be patrolling the area which that route represents. Figure 3 displays a general idea of 


smuggling routes as interpreted by the Drug Enforcement Agency. 
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Figure 3 Drug flows from South America into the United States. 
Source: DEA Map | - 1987 (Revised 1989) 


22 


A third assumption is made when computing the probability of detection on 
the routes. Each asset represents a sensor which adds a probability of interdiction to the 
route that it covers. It is also assumed that an asset covering a leg of the route acts 
independently of other assets on that and other legs, and that detection on any leg of a 
route can be taken as an event independent of all other legs on that route. Realistically, 
a patrol aircraft on a route, which detects a suspicious surface contact 1s likely to be in 
either direct or indirect contact with a surface asset on that leg or a subsequent leg, giving 
information about the presence of that contact to the surface asset and increasing the 
surface asset’s probability of detecting and interdicting that contact. Similarly, a surface 
asset may detect a suspicious aircraft and report it, etc. For modeling purposes, however, 
we believe these assumptions to be reasonable. Independence between assets and routes 
will result in SOAR providing an upper bound on the measure of effectiveness. This 
consequence should be considered when estimating the probabilities of interdiction input 


to the model for each asset type against each method of smuggling. 
3. Data Structures and Computations. 


a. Smuggling Routes. 
The representation of routes is accomplished using a_ hierarchical 
adjacency list of geographical points. The geographic points are listed in a data file 
similar to that shown in Appendix A Section A.5., which includes points of origin for 


drug shipments, intermediate route points, and points of entry into the United States. 
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The information for the routes themselves is taken trom another data file. 
A description of the inputs and an example of the data file are in Appendix A Section 


A.3. 


b. Assets. 

Appendix A Section A.4. describes the input information and a data file 
for interdiction assets. Although we have listed assets as SHIP], PAIR], etc., asset need 
not necessarily refer to a single ship, aircraft or ground unit. An asset may be considered 
whatever mix of elements 1s required to achieve the range and probability of interdiction 
input to the model. For simplicity, we will continue to refer to assets as single units. 
Asset types are used to describe assets with the same range and probability of interdiction. 
By assigning additional asset types, the matrices containing this information could be 


expanded to include unique input data for particular assets. 


c. Computation of route probabilities. 

A shipment of drugs will be interdicted on a particular route with a 
probability computed in the subroutine PROBCOMP. In this subroutine, the legs of each 
route are checked to determine if assets are located such that the leg falls within the range 
of the assets. Up to this range, the asset contributes a constant probability of interdiction 
to the route. Beyond this range, the probability decays exponentially as a function of the 
additional distance. The probability of interdiction on a route is then computed as simply 


the combination of independent probabilities for all assets covering all legs of the route. 
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The probability decay function chosen results in probabilities for 
individual assets decreasing to approximately 13.5% of the onginal value at twice the 
assetrange. For example, an asset with a given constant probability of 0.10 to 100 miles 
will have a probability of interdiction equal to 0.0135 at 200 miles. This method of 
computing probabilities can be supported by considering the given range of an asset as 
that area which can be patrolled with a constant probability, and exponential decay of the 
probability representing the decreasing effectiveness of the patrolling as the distance 


between the asset and the route increases past that range. 
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IV. THE ANNEALING MODEL 


A. GENERAL DESCRIPTION 

The simulated annealing method is so-called for its resemblance to the annealing 
of solids, in which a material is heated to the point where all particles arrange themselves 
randomly in the liquid phase, and is slowly cooled so that all particles arrange themselves 
in the low energy ground state of a corresponding lattice. [Ref. 22: p. 2] 

The method starts at either a predetermined or random initial point. In the feasible 
space of the independent variables, a random walk samples the objective function 
[Ref. 23: p. 210] (analogous to the energy stale of a _ physical system). 
Each step providing an improved objective function value is automatically accepted. 
Steps which are detrimental to the objective function value are accepted according to a 
Boltzmann probability tunction, p = exp(-AE/K,T), which is dependent on the energy and 
temperature of the physical system. This allows the path to walk out of local and global 
optimal points [Ref. 23: p. 211]. The process is often called a Metropolis loop, named 
for the author who used a model simulating the evolution of a solid to thermal 
equilibrium in 1953 [Ref. 24: p. 14]. Obviously, the algorithm could run 
indefinitely, walking in and out of the global optima. One way to prevent this from 
happening is to specify the number of steps the algorithm is to take. However, this in no 
way guarantees reaching the global optimum. Recently much research has been put into 


developing cooling or annealing schedules. Essentially, such schedules decrease the 
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probability of accepting a detrimental step by decreasing the temperature of the system 
as the global optimum point is approached. The random walk is terminated when a 
predetermined number of successive tals does not produce an acceptable objective 


function value. [Ref. 23: p. 212] 


B. APPLICATIONS 

The simulated annealing method has been applied to optimization problems in 
numerous fields, including physics, operations research, numerical analysis, biology, 
materials science, game theory, code design, etc. [Ref. 22: p. 14]. 

There are many advantages and disadvantages associated with the simulated 
annealing method. The algorithm is applicable to many different opumization problems 
and is easy to implement. The greatest disadvantage is the potentially prohibitive amount 
of time required to converge to a near-optimal solution. This is dependent on the cooling 


schedule and the step size used in the random walk. [Ref. 25: p. 91] 


C. THE GENERALIZED ANNEALING ALGORITHM 

As stated previously, simulated annealing algorithms have been used almost 
exclusively in applications which minimize the objective function. The probability 
function is especially well-suited to cases where the objective function minimizes to zero, 
or where the optimal value is known. In these cases the probability of accepting a non- 
improving step is driven to zero as the optimum is approached. 

In the case where the optimal value is not known, the typical algorithm calls for 


selecting some initial estimate of the optimal value (@,,) and allow the random walk to 


2/ 


sample the objective function value (), continuing until o - o ,, becomes negative, then 
decrease ,, repeating this process as required. On the other hand, if the initial estimate 
of ,, 18 too low, it is increased as part of the cooling schedule, so that  - o,, will 
approach zero. [Ref. 23: p. 213] 

We have attempted one of many possible applications of this method to a 
maximization problem, by taking an initial guess for the optimum as an input value, as 
well as the number of steps to be taken before reducing the guess. Cooling the guess 
(0,,) 18 accomplished by keeping track of the highest value found so tar (@,), and 
reducing @,, according to the function 6,,, = 0, - @ * (0, - 0s), where @ is a 
predetermined fraction which controls the cooling speed. This ensures that 0,, is never 
reduced past the best value found in the random walk. If the random walk produces a 
value which is greater than either the current 6,, or the input guess, then both are 
increased a significant amount above the higher value (we have chosen @,, = 1.5 * 5). 
The cooling of @,, again proceeds according to the previously descnbed method. 

In our application of the algorithm, we begin with a large step size (Ar) in order to 
sample a significant portion of the feasible region. As 0,, 1s reduced (or cooled) we also 
reduce the step size to a minimum value (Ar,,,,) by a similar equation, sampling a tghter 
area of the region. When a potential point of convergence is found, Ar is reduced to a 
small value (Ar,,,) in an attempt to find an acceptable move in the very close vicinity. 
If an acceptable move is found, the step size is reset to the previous value (Ar,,,), and the 


annealing continues. 
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This algorithm is not guaranteed to locate the objective function value associated 
with a global optimum. The value 8 in the probability function controls the rejection 
ratio, and will also determine, in part, how close to the optimal value the algorithm 
converges. For very low values of B, the algorithm will converge at, or very close to, the 
optimum, when the optimum value is known. 

The feasible space and step generation used for our model are very simple, 
consisting of upper and lower limits on the latitude and longitude for asset locations, and 
random steps with size Ar. For example, we did not attempt to prevent assets from 
moving onto land, restricted airspace, or within search range of other assets. These are 
the more detailed geographic, operational, and political constraints, such as those 
discussed in Chapter I, Section B, which could be used to control possible asset moves. 
Such constraints could be represented in a subroutine that tests feasible locations and a 
set of rules controlling the possible moves of assets relative to each other. 

Our algorithm, shown in Figure 4, is very similar to that presented by Bohachevsky, 
Johnson and Stein for a minimization problem allocating ballistic missile interceptors 
(Ref. 26]. The following terms used in the algorithm are defined: 


¢ @ is the objective function value. A subscript of 0 indicates the current value, | 
indicates the value at the new location. Other subscripts are as defined above. 


¢ ) is the feasible space, that is, the set of latitudes and longitudes where assets may 
be located. 


¢ x is the set of latitudes and longitudes for all mobile assets. A subscript of 0 
indicates the current location, | indicates the new location. 


¢ Ar is the step size. Subscripts min and low indicate the minimum step size used 
when cooling and the step size desired in attempt to find a more precise solution. 
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Step 1: Compute the objective function value for initial locations, 0p. 
Set $,, = input guess. 


Step 2: For each mobile asset, generate two independent standard normal random variates, 
A, and A,, and compute the components of the random direction vector d: 
d= A/(A,**2 + A,**2)"*0.5, i = 1,2. 


Step 3: Compute candidate locations for assets, set x, = X, + Ar * d. 
Step 4: If the candidate locations are infeasible, ie, x ¢ 22, repeat Step 2. 
Step 5: Compute the new objective function value, >,, and set Ad = , - %o. 
Step 6: If the step increases the objective function value (Ad 2 0), accept: 
set k= 0, xX, =X,, 0, =6,, andn=n+1. 
If the step size has been reduced in Step 8 (k > S), reset step size Ar = Ar,,,. 
lf >, > 5, Set O, = 0). 
6a: Increase 9, and initial guess if new objective function value is greater, ie, 
lf 0 > 0 scel@ = — shoe 
If d, > initial guess, set guess = 1.5 * 9,. 


Set n = 0. 


6b: If N accepted steps taken (n > N), cool o, and the step size: 


set ,, = $,,- a * (0,, - $,) and Ar = Ar - a * (Ar - Ar,,,). 
Set n = 0. 
Go to Step 2. 


Step 7: If the step is a decrease (Ad < 0), compute probability of acceptance and 
determine whether to accept or reject. 
Set p = exp( B * Ad /(o,,-,)). Generate a uniform random variate V. 
If V > p, reject, set k =k + 1. If V <p, accept, set k = 0, X = X,, 0) = 0). 


Step 8: If the stopping rule has not been reached (k < S), go to Step 2, 
otherwise, if k = S, attempt to find a more precise solution with small step size, 
save current step size, Ar,,, = Ar, set Ar = Ar,,, and go to Step 2. 

Step 9: If the final stopping rule has not been reached (k < 2 * S), go to Step 2, 


otherwise stop and return $, as the optimal solution found. 


Figure 4 Annealing algorithm to optimize interdiction asset locations relative to a 
stochastic network. 
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¢ is the control parameter for the probability function (analogous to Boltzmann’s 
constant). 


¢ ais the fraction for reducing o,, and Ar 
¢ N is the number of accepted steps at each value for @,, 


¢ S is the stopping rule, the control for the number of attempted moves rejected 
before the solution is declared found. 


¢ k is used to count the number of consecutive rejected moves 
* nis used to count the number of accepted moves for each level of 9,, 


The FORTRAN code for the model can be found in Appendix B. 


D. RESULTS 

The optimization model developed in this thesis was applied to a smuggling 
network of eleven air and sea routes in the Caribbean region. Figure 5 displays this 
network, formulated from the sample data sets in Appendix A. Some areas of the routes 
used are intuitive, for example, those running through the Yucatan Channel and the Mona 
and Windward Passages. However, these are meant to be illustrative of actual routes, and 
are not based on intelligence or other official information. The model was run on a 
number of different computers, SUN SPARC 1+ and 2 workstations, Tektronix 4336, 
Silicon Graphics IRIX, and an AMDAHL 5990-500 mainframe computer. The four 
UNIX systems yielded identical results, while the VMS-based mainframe computer found 
similar results to the UNIX systems. 

Two applications of the model were tested. The first examined the locations at 


which the annealing algorithm converged. Ten runs each were performed with five 
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Figure 5 Sample Caribbean smuggling network 


mobile assets, three ships and two aircraft. The resulting locations from one run are 
plotted in Figure 6. It 1s clear that some of the locations could easily be improved. This 
is due in part to the value of the probability function control parameter, B, used to ensure 
convergence of the algorithm in a reasonable amount of time. Each of the ten runs found 
different final locations, but all appeared reasonable and the optimal values found were 
similar, with mean equal to 75153.2 (dollars x 1000) and a standard deviation of 7204.1. 

The second application explored the response of the model to a changing number 
of assets. Each run used equal numbers of ships and aircraft. Parameters of assets are 


from the sample data set in Appendix A. We found that the function, 


a2 


ASSETS 
@ = Patrol air 


@ = Ships 


Three ships, two aircraft 
Optimal value = 78047 





Figure 6 Results of five asset run (one of ten) 

d = 615430.6 - 71.07 * (80 - ASSETS) ** 2.07, determined using non-linear regression, 
fit the data well (R* = 0.9922). This particular function is dependent on the number of 
assets at the last data point, and 1s extremely simplistic, with total numbers of assets as 
the independent variable. Runs with more than 80 assets will change the parameters of 
this function, but will not significantly increase the maximum. Although the fit of the 
curve is good, the variance at each level of assets cannot be assumed equal, mainly due 


to the small sample size at the levels (2). 


The regression model itself is not a fundamental result, but can be used to obtain 
useful information. With 80 assets we are very nearly saturating the routes, obtaining 


little additional benefit from adding assets. We will assume for the following discussion 
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that 80 assets induce the highest possible total cost to the smugglers. In Figure 7 we have 
plotted a curve depicting the fraction of maximum value achieved versus assets. This 
fraction represents the regression function value divided by 615430.6, the function 
evaluated at 80 assets. From this, we see that 50% of the maximum attainable effect on 
the smugglers is obtained with 24 assets, 75% with 36 assets. In other words, 75% of the 
maximum effectiveness 1s achieved with fewer than half the maximum number of assets. 
It is highly unlikely that such quantities of assets would be available for employment at 
any particular ume. However, these results are entirely dependent on the specific 


parameters used and our simple asset mix. 
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Figure 7 Fraction of maximum value found (at 80 assets) vs assets 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

In this thesis we have demonstrated that simulated annealing can be employed to 
find near-optimal locations for interdiction assets relative to a network with stochastic 
flows that are determined by a simulation. Our approach also demonstrated an application 
of this type of model to explore the relationship between the capabilities of interdiction 
assets, their number, and the resulting effect on the overall effectiveness of the 
interdiction effort. Although our model uses a simple representation of interdiction assets, 
combinations of assets, and the geographic, operational and political constraints which 
would normally be used to restrict feasible asset locations in the drug interdiction 
problem, we believe that the results demonstrate this methodology’s viability for 
development into the type of analytic tool described in Chapter I, Section B. 

The primary obstacle to further use of this methodology in counter-narcotics 
operational analysis and planning is the questionable reliability and limited amount of 
information available on smuggling activities. The credibility of results from this type 
of model, including the validity of SOAR, is highly dependent on the accuracy of 
information used in constructing the route network, and assigning costs to smugglers and 
parameters to assets. The suspect reliability of smuggling data will make it difficult to 
assess this credibility, short of actually employing operational assets and observing the 


resulting effects on smuggling and drug use. Even this level of verification would not 
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necessarily provide an accurate measure of credibility since the targets of drug interdiction 
are the only true source of the information required, and they are not likely to come 


forward with such details. 


B. AREAS FOR FURTHER RESEARCH 

The next logical step in the development of an analytic tool useful to the counter- 
narcolics program is to apply the methodology explored in this thesis, or the model 
presented itself, to an expanded problem. Such an expansion should include all primary 
cocaine smuggling regions (Atlantic/Caribbean, Pacific, and Southwest Border), 
constraints reflecting the geographic and political considerations for asset locations, and 
Operational inputs, such as budget feasible asset combinations and capability restrictions 
of assets. These constraints could easily be included in a subroutine which tests proposed 
asset moves for feasibility. Analysis of the results from an expanded problem will 
provide a more realistic evaluation of the utility that may be derived from this type of 
tool. 

Expanding the size of the problem produces a major drawback in the use of the 
model developed in this thesis, extensive computer run times. With a network of only 
eleven routes, the five asset runs described in our results averaged seven hours each on 
the SUN SPARC 1+ workstation. Adding routes to the problem will significantly 
increase the run time of SOAR. Furthermore, our testing used only one trial of SOAR 
to obtain each new objective function value, while future testing should use the average 


of multiple trials to reduce the variability of the simulation results. This will increase run 


times even more, by approximately a multiple of the number of trials. For these reasons, 
analysis of the SOAR model to optimize performance, or translation into a more efficient 
programming language, and the use of high speed computers would be beneficial to future 
efforts. 

A final consideration involves the variety of measures of effectiveness (MOEs) 
being used by drug interdiction agencies, as discussed in Chapter III, Section B.1. 
Development of a well-defined, quantifiable MOE for the entire drug interdiction program 
would be very valuable to future analysis. The MOE produced by the SOAR model for 
Our Optimization, total cost to the smugglers, relies heavily on information which may not 
be particularly accurate, as would be the case for any MOE that reflects the goals of the 
interdiction program as a whole. Continued intelligence efforts to gather accurate 
information on illegal drug smuggling activities and improvements in the timeliness of 
such data will help in defining an overall MOE and assessing the credibility of results 
from analytic tools used in counter-narcotics planning, and will make the drug interdiction 


program more effective. 
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A. MODEL INPUT AND OUTPUT 
INPUT DATA SETS 
1, ANNEAL DATA 


a. Variables 


Variable name Description 


STEPSIZE Step size for the annealing 
algorithm in miles. Real. 


STEPSIZEMIN Minimum Step size for the 
algorithm in miles. Real. 


STEPSIZELOW Step size in fine search 
for optimum value in miles. 
Real. 


BETA Control parameter for 
annealing. Adjusted to 
control the 
acceptance/rejection ratio. 
Real. 


S The stopping rule for the 
algorithm. This number of 
consecutive rejections stops 
the algorithm. Integer. 


ALPHA Parameter used to decrement 
the estimate of the optimal 
value and step size. 

Adjusted to control the speed 
of cooling. Typically between 
0.01 and 0.2. Real. 
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GUESS 


BIGN 


b. Sample data set 


Inital estimate of the 
opumal objective function 
value. Real. 


The number of steps to be 
taken at each level of 
cooling. Integer. 


INPUT PARAMETERS FOR THE ANNEALING ALGORITHM 


SlEPSIZE (MILES) 
Pea! STEP SIZE 
ewe STEP SIZE 
CONTROL PARAMETER 
SlOPPING RULE 
ALPHA 

GUESS AT MAXIMUM 
ere PS PER LEVEL 


100000000.00 
1000 


2. SOAR DATA 


a. Variables 


Variable name 


NUMMETHODS 


NUMPHASES 


NUMROUTES 


NUMTRIALS 


ENDTIME 
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Description 


Number of methods for 
smuggling drugs. Integer. 


Number ot phases. Integer. 


Number of routes available. 
Integer. 


Number of times the analysis 
period is to be simulated. 
Reported results will be 
averaged over the trials. 
Integer. 


Number of days to be 
analyzed, not including the 
runin period.Integer. 


RUNIN 


[X 


DRUGCOST 


EXSHIPMENTINTVL 


EXSHIPMENTSIZE 


METHODNAME(M) 


CAPCOST(M) 


CAPACITY (M) 


Number of days to run the 
simulation before commencing 
the analysis. Integer. 


Seed for random number 
generator. Integer. 


The cost at the source of a 
kilogram of cocaine. Real. 


Expected time between 
shipments. Real. 


Expected shipment size. Real. 


The name of drug smuggling 
method M. Character. 


The cost to the smuggler of 
an unsuccessful shipment 
using method M. Does not 
include the cost of the drug 
or the cost associated with 
the route. Real. 


Maximum amount of drug that 
may be shipped by method M. 


Real. 
b. Sample data set 

METHODS PHASES ROUTES TRIALS RUNIN ENDT IME MEMORY 

3 2 2 rt) 60 21 
DRUGNAME DRUGCOST SHIPMEAN SHIPSIZE DAILY AMT 
COCAINE 50050 Carel 250.0 S500 
METHOD RISK COMP RISK EXP SEIZE COST MAX SHIP 
AIR 1200000.0 2.0 200000.0 2000.0 
SEA 1600000.0 2.0 40000.0 16000.0 
LAND 10000.0 220 5000.0 50.0 

01 

61 
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3. ROUTES DATA 


a. Variables. 


Variable name 


ROUTENAME(R) 


ROUTECOST(R) 


ROUTEMETHOD(R) 


ROUTELOC(R) 


FIRST(R) 


NEXT(R,P) 


Description 


The name of route R. 
Character. 


The cost to the smuggler for 
using route R, incurred 
whether or not the shipment 
is successful. Real. 


The index of the smuggling 
method that is used on 
route Ry intever 

(l= Air, 2=Sea-3=Lana): 


The index of the geographic 
area of route R. Integer. 
(1=Atlantic,2=Pacific, 
3=Southwest Border). 


The index of the first point 
of route R. Integer. 


The index of the point 
following P on route R. 
Integer. 


b. Sample data set. 


ROUTENR ROUTENAME ROUT ECOST ROUTEMETHOD ROUTEBEGe 


POINTS 

1 MONA PASS 160CU20u 2 1 
1 10 ial 38 

2 WINDWARD 16000.00 2 il 
2 12 13 14 Sy) 

3 YUCPLAI 16000 .00 2 il 
fe ie Sie 

4 YUCALA] 16000 .00 gy 1 
S le a2 

2 YUCTE Pe OU e010 2 1 
5 ilies) oe 

6 YUCFLA2 16000 .00 2 1 
4 is) 36 

7 YUGEIEAS 16000 .00 2 1 
4 ib. 35 

8 YU Ez TSEViOCz00 2 I. 
4 15 og 

) CUBA 20000.00 1 i! 
3 Shs 

10 YUCALA2 20000.00 | 1 
7 34 

ml YUCTEX3 20000.00 1 ] 
a sill 
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4. ASSET DATA 


a. Variables. 


Variable name 


NUMASSETS 


NUMASSETTYPES 


ASSETNAME(A) 


LATASSET(A) 


LONASSET(A) 


Boor TTY PEA) 


ASSETLOC(A) 


ASSETRANGE(AT,M) 


IPD(AT,M) 
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Description 


Total number of assets, 
mobile and fixed. Integer. 


Total number of asset types. 
Integer. 


The name of asset A. 
Character. 


The latitude of asset A. 
Real. 


The longitude of asset A. 
Real. 


Index of the type of asset A, 
for example, 

| =air,2=ship,3=land, 
4=aerostat balloon. Integer. 


Index of the geographic 
location of asset A. Same as 
ROUTELOC. Integer. 


Range for asset type AT 

versus smuggling method M for 
which the probability of 
interdiction 1s a constant. 

Real. 


Probability of interdiction 

for asset type AT versus 
smuggling method M on a 
route, which is constant out 

to range ASSETRANGE. Real. 


b. Sample data set. 


TOTAL NUMBER OF ASSETS = 8 
NUMBER OF ASSET TYPES = 4, PATROL AIR,SHIPS,LAND UNITS, AEROSTATS 
ASSET NAME LATITUDE LONGITUDE TYPE LOCALE 
1 PAIRI 16.00 75 85 1 1 
2 PATR2 17.45 80.95 1 1 
3 SHIP1 19.00 68.00 2 1 
4 SHIP2 19.00 75.00 2 1 
5 SHIP3 21.98 86.18 2 1 
6  CUDECH EEL 247 4Ad 81.30 4 1 
y MARFA TX 30.19 Oecd 4 3 
8 GREAT EXUMA Denes 54g 4 1 
ASSET RANGE VS METHOD 
METHOD 
ASSET AIR SEA LAND 
AIR 100.0 100.0 25.0 
SHIPS 50.0 50.0 Dee 
LAND pope 0.0 BAG 
BALLOONS 200.0 0.0 0.0 
P(INTERDICT) ASSET VS METHOD 
METHOD 
ASSET AIR SEA LAND 
AIR 0.1 al or) 
SHIPS Ome On 0.0 
LAND 0.1 0.0 Osa 
BALLOONS C8. 0.0 0.0 
5. GEOGRAPHIC DATA 
a. Variables. 
Vanable name Description 
POINTNAME(P) The name of geographic 
point P. Character. 
LAT(P) The latitude of point P. 
Real. 
LON(P) The longitude of point P. 
Real. 
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b. Sample data set. 


POINT POINTNAME Ey OI B) LONGITUDE 
1 PT GALLINAS Zs (ele Se. 
2 RIO HACHA bipalre S (0, oO 
3 SANTA MARTA ST heels: TA 10 
4 TURBO S200 76.45 
5 BUENAVENTURA 28 das 0 
6 TUMACO 45 TE.45 
7 MEDELLIN 6.10 (eeee se 
8 BOGOTA a aes 
2, CALI S2aZ0 T6328 

10 SMONAPASS NE re C6700 
ial NMONAPASS 20,50 68.00 
2 SW INDWARD ices: ewe l le! 
IIgs NW INDWARD 20330 Toee7 
14 SI RALCH IS Oe Be: Focaech| 
Is YUCATAN OA ees) eo crrae Les: 
iis CENT MEXICO 25.00 300 
ig NOGOLES AZ Shee 4 Ns oe 
18 BISBEE AZ Ses) Lo 5 
ili, DOUGLAS AZ SD Es 0 
Z0 YSLETA TA Se, 106.18 
21 PRESIDIO. IA Zo Gas 2 3 
ae DELRITO.TX 29%.2 1 Ors: 2 
eS DOLORES TX Zea? 9934] 
24 BROWNSVILLE TX ZocoS hi ieee 0 
as) SAN DIEGO CA e241 TRGe 57 
26 CALEXICO CA e241 Le 0 
20 Fara NM ee a8, Oe. 5 
28 ee2 NM ante OS neo 
a9 PREEPORT atx Clots), owe | 
30 GALVESTON TX ZOeLS UA a6 
aril HOUSTON TX 25 6 oo ea 
2 PORT LAVACA TX Z3a36 10538 
a3 SABINE TX 29.44 Gee 4A 
34 MOBILE AL 20542 86.03 
3s) PENSACOLA FL a0 a25 oie oral le 
e6 Serer is Coals ore ae: 
of FLAMINGO FL 2220 80.255 
38 DAYTONA BCH FL 7s aad Bh ee 2 
oo MIAMI FL 25.45 odes! 
40 LOS ANGELES CA 34.00 ies ed as, 
4l BAJA TIP 22.64 OS 82 
LD: S MEXICO ES Ws We: cs melles! 
43 Pape Tet L3432 1052315 
44 PAGTiELCZ 20200 11S 200 
45 PAG TRIG 3 20 a Osa) 
46 PACIPIC4 Loo POS. U0 
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B. 


OUTPUT DATA SETS 


I. SOAR OUTPUT 


Output of the model includes an echo of the input parameters, and the 


following Statistics: 


L. 


Y 


2 


Expected attempts per trial. The average number of shipments attempted during the 
period being analyzed. 


Expected successes per trial. The average number of successful shipments during the 
analysis period. 


Expected interdictions per trial. The average number of unsuccessful shipments 
during the analysis period. Expected attempts = Expected successes + Expected 
interdictions. 

Success rate. The proportion of shipments that were successful. 

Interdiction rate. The proportion of shipments that were unsuccessful. 

Cost of incomplete shipments (in thousands). The average cost to the smuggler 
because of unsuccessful shipments, including the cost of the method (such as an 
airplane), the cost of the drug, and the cost of the route (such as gasoline). 

Cost of completed shipments (in thousands). The average cost to the smuggler 
because of successful shipments, including the cost of the drug and the cost of the 


route: 


Total cost to smugglers (in thousands). The sum of the cost of incomplete shipments 
and the cost of completed shipments. 


For each drug, the quantity that the smugglers attempted to ship. 


10. For each drug, the quantity that the smugglers successfully shipped. 


11. For each drug, the quantity that the smugglers lost because of unsuccessful shipments. 


12. For each route, the expected attempts, successes, and failures are reported. 


13. For each phase and each route, the expected attempts, successes and failures are 


reported. 
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14. For each drug and each route, the expected quantities shipped, captured, etc.[Ref. 2: 
pp. 42,43] 


2. ANNEALING MODEL OUTPUT 


ies TOTAL COST TO SMUGGLERS = 78047 .95 
(IN THOUSANDS ) 


Orl VALUE FOUND IN o63 UNSTERS 
SS i@o EPS NOT ACCEPTED 


deo INPEASIBLE MOVES FOUND 


FINAL LOCATIONS 
ASSET LATITUDE LONGITUDE 


PATR1 Pee ans) Garo 4 
PATR2 Zoo ieee oe: 
Sie ea PH hs oe eis 
Sie Ono oS 86.613 
SiGes TA Sia 6077279 
SUDLOE KEY FL 24% 4 100 SPes00 
MARFA TX Bi) 2190 10400 
GREAT EXUMA 23.350 J3-4 70 
foes PROBABILITY OF INIERDICTION BY ROUTE 
ROUTE ROUTENAME PROGINTERDICT 
i MONA PASS Oi 7 4 
Z WINDWARD OZ 130 
3 YUCFLA! Oars FAO KS 
4 YUCALA1 On 2176 
5 OC Ee I O23 146 
6 VUGrEAZ O22) 
7 YUCE LAS One 9 
8 TOC TE 2 Vn 2519 
9 CUBA ee 249 
10 YUCALA2 G22909 
ial YUGTE Ss On F692 
Sey TIME = 17 704-3 SECONDS 
MAX VALUE FOUND = Toole EeO07 
Reni GUESS VALUE = S710 84SE+07 
Nba LOC SEED = 295007429 


mine COMPLETED 
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* 


* 


* 


4 + 


+ + 


+ + 


APPENDIX B. FORTRAN CODE 


PROGRAM ANNEAL 


WRITTEN BY LP JAMES +d. HENRY i asl) 
NAVAL POSTGRADUATE SCHOOL, CODER oc. 
OPERATIONS ANALYSIS CURRICULUM 


PROGRAM VARIABLES 


INCLUDE ‘COMMON .FOR’ 
REAL PROBMOVE, 
TLAT (MAXASSETS) , TLON (MAXASSETS), 
OBJVALO, OBJVAL1, OBJVALM, OBJVALB, DELOBJ, 
VP TOrCOs! , CPUrinig. 
A(2) DENOM TARR A ei. 
INTEGER K, bLITILEN, 1, ROUTENK, MOVES ,NCIMOVE hinEsAs 
LOGICAL FEASIBLE, ACCEPT, PINTS soorinirog 


Sg a Sig i 


CALL SUBROUTINE TO READ IN THE INITAL DATA 


CALL GETDATA (GOODINPUT) 
ISA d es, Di I. 
bSBED lS 


TE. NOT 2 GOocbin rule ler 
OBJVALM = GUESS 
DO 50 l=1,NUMASSETS 
TELAT (1) =) Garesser ei} 
TEON(..)> = ShONRSSEr ap 
50. CONTINUE 
RUN THE SIMULATION WITH INITIAL Daly eortiis 
CALL SMUGSIM( TLAT, (LON TOTeccer 
OBI VAnOw= “TOrecosl 


OBJVALB = OBJVALO 
1F (OBJVALB.GT. "GUESS )@ CUbS soe eel) ee 


BEGIN THE ANNEALING 


MOVES = QO 
NOTMOVE = 90 
INFEAS = 0 
PINISH-= ~FPALSE. 
Ke=-0 


od TGs 0 
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POO CCON TINUE 


PINS tw LOCATIONSEOR ALLuASSETS BACEPT 
BALLOONS (TYPE 4) 
DO, 140 Is] WeNUMASSETS 
LEC ASSEV IE) NE. 4) THEN 
rz 0 CONT INUE 

CALEVUNCKPCUISEEDI,, A,e2} 
DENOM == “SORT (AN 1) **2 455A ( 2) **2)} 
THAT (1) = UATASSET (i) + (DELR*A{1.)./DENOM 
TLON(I) = LONASSET(I) + DELR*A(2)/DENOM 


CHECK PORK FEAST BIUILy 

PE ASSET EOC 1)7. EO? 1.) THEN 

CALE LANTCHECK{(TLAT (1), TLON (1), FEASIBLE) 
ELSEIPTASSETLOC(L): 50. 2) THEN 

CABL PACCHECK(TLAT( 1); TOON). PEASIBLE) 
ELSE 

CALL SWECHECK (TLAT (1) TEONV) PEASIBLE} 
ENDLE 
LEC NOTTS PEAS EBLE MINEEAS?— UNE EAS ee 


LE JNOT. SFEASIBLE) GO. TOR iTZ0 
ENDIF 
140 CONTINUE 


RUN THE SIMULATION WITH NEW DATA POINTS 
CALL SMUGSIM(TLAT, TLON, TOTCOST } 
OBOY Alii = TOreCOsT 
DELOBI = OBJVAII — OBIVALO 
ADJUST MAX VALUE GUESS IF REQUIRED 
IF{(OBJVALI .GT. OBJVALB) OBJVALB = OBJVALI 
IF{OBIJVALB «Gl. GUESS) GUESS = 1.5 * OBJVALB 
IF(OBJVALB .GT. OBJVALM) OBJVALM = GUESS 


DETERMINE WHETHER OR NOT TO ACCEPT THE MOVE 
Tee lWOBIS GE. "6205 THEN 
ACGERT =~. TRUE. 
Boo 
CK = OBJVALM - OBJVALO 
Le Ck eo a0. Oy Sek a=. 000007 
Pea EAP Cee ra. * UBLOB) <7 CK) 
CALL ERNDPGYISEEDL, V1) 
[iP PGE. Vy) wen 
ACCEL = LeU E.. 
PlosSk 
ACCEPT = .FALSE. 
K seas 1 
NOTMOVE = NOTMOVE + 1 
ENDIF 
ENDIF 
te (Ae CE rr). THEN 
COUNT = 0 
DO 160 70-= 1 NUMASSETS 
GATASSET (1) = TLAT( LE) 
LONASSET(I) = TLON(I) 
160 CONTINUE 
MOVES = MOVES + 1 
OBUVALO -="OBJVALI1 
Pritben = LIPPoEN + 1 
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ENDIF 
10 le OR Se laa 
DELROLD = DELR 
DELR = DELRLOW 
ELSEIP(K . CGE a2. 2S mee 
FINI St e=— 4 RUE 
ELSEIP(K .GT. -S .AND SSAC Eee wie 
DELE -=]DELROLD 
ENDIF 


IP ( LbITToRN GE. BIGh) are 
OBJVALM = OBJVALM - ALPHA * (OBJVALM - OBJVALB) 
DELR. = DELKR-= AGPHA > (DEER — sien eh) 

ENDIF 

IF (OBJVALM .LT. OBJVALB) ObdVALM}="CBIVALB 


PF (.NOT.SFINTSH) 2 GO, Te wu 
CLOSE S5) 


WRITE THE RESULTS 
OPEN( 75, FIbe— ANNEAL AOUI) 


WRITE (75,520) OBUY FEO Ue) 

520 FORMAT (1X, “MAX TOTAL COs) 10. SMUGGLERS. =n iis 
WREVE( 75.530) 

530 FORMAT ( 10%," (IN THOUSANDS =e 


WRITE(¢ 3S, sUU) MOVES 
500 FORMAT (IX, “OPT VALUE FOUND MIN 7a Slee | 


WRITE (73, 501} ieiMoveE 
SOT PORMAT( 19%, LIU s Sale Rovio ea eer ars 


WRIETE( 7S, 73 
WRITE 75, 00sUNE EAS 

910 FORMAT (IX, 1105- INFEASIBEEE, Mo Peer eunh 
Web EA 


WRITE oe ? 
WROTE) S57) 
WRB 5954) 

034 FORMAT( 1124," FINAL LOCATIONS. 

WRIStE 7 5,535) 
535 FORMAT(1IX, “ASSET. bDATITUDE 77 heer iVbias) 
DO So0ylai_ »NUMASSETS 
WRITE(75,540) ASSETNAME (I) , LATASSET (I) , LONASSET (T) 
540 RORMAT (yeaa W054 ba Ome 
550 CONTINUE 


WREPEM oy) 

WRI 755 eo) 
589 FORMAT(’FINAL PROBABILITY OF INTERDICTIONGE aa Our.) 

WRITE (7 55.52 0)) 
590 FORMAT {( ‘ROUTE’ ,5X, 'ROUTENAME 74x) PRep ner piel) 

DO 600 ROUTENR = 1, NUMROUTES 

WRITE (75,610) ROUTENR, ROUTENAME (ROUTENE = 

& PROBINTERDICT {ROUTENEG ss) 
610 FORMAT (1%, 14, 5%, Al07 5x7) oeeD 
600 CONTINUE 
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Se a = DTIME(TARRAY } 
WRITE Cr, * ) 
WRITE( TS, f)  SeRU iii ess | CPULIME,.” SECONDS” 
WRITE(75,*)‘ MAX VALUE FOUND = ‘,MAXNOW 
WRITE (7 Se)! FANAL GUESS VALUE = * ,OBJVALM 
eae *) 
Vivian) NE CP L@GeShED = .*) ISEEDI 
WEE, 7) RUN COMPLETED” 
CIGOSE (75) 
SLOP 
END 

aa. SUBROUTINES TO CHECK FEASIBLE REGIONS ee, rn nn eee 


SUBROUTINE LANTCHECK (LAT, LON, FEASIBLE) 


Reais bAT GON 
LOGICAL FEASIBLE 


PeAoleGk = SPARSE. 

i GON bs. Ofe0e.AND2) LON. .GE267.0 SAND: 
& Ate by . 2omdes. AND. UAT Ghee 142.0) PRAGIBBRE. = ~TRUE.. 
RETURN 

END 


KEKE KKK KEKE KEKE KKK KKK KE KKK KEK KKK KEKE KKK KKK KKK KEKE KK KKK KE KEK KEKE KK KKK KK KKK KKK KKK 


SUBROUTINE SPACCHECK LAT LON, FEASIBLE) 


REAL LAT, LON 
BOGICAL FPEASIBLE 


FEASIBLE = .FALSE. 
* 
Pei LATeebE. 0.0) THEN 
FEASIBLE = .FALSE. 
ELSEIF(LON .LE. 120.0 .AND. LON .GE. 118.691) THEN 
IF(LAT .LE. 33.773) FEASIBLE = .TRUE. 
ELSEIF(LON .GE. 118.0) THEN 
IF(LAT .LE. 32.50) FEASIBLE = .TRUE. 
ELSEIF(LON .GE. 115.50) THEN 
IF(LAT .LE. (-132.70 = 1.40*LON)) FEASIBLE = . TRUE. 
ELSEIF(LON .GE. 113.00) THEN 
IF(LAT .LE. (-64.40 + .80*LON)) FEASIBLE = .TRUE. 
ELSEIF(LON .GE. 107.00) THEN 
IF(LAT .LE. (-69.1667 + .8333*LON)) FEASIBLE = .TRUE. 
ELSEIF(LON .GE. 84.50) THEN 
IF(LAT .LE. (-39.4444 + .5556*LON)) FEASIBLE = .TRUE. 
ELSEIF(LON .GE. 79.00) THEN 
IF(LAT .LE. (-15.5455 + .2727*LON) .AND. 
& LAT .GE. (121.5 - 1.5*LON)) FEASIBLE = .TRUE. 
ENDIF 
RETURN 
END 


KKK KKK KKK KKK KEK KKK KKK KKK KKK KKK KKK KKK KEKE K KKK KK KKK KEKE KKK KKK KKK KEK KKK 


SUBROUTINE SWBCHECK (LAT, LON, FEASIBLE) 


REAL LAT, LON 
LOGICAL FEASIBLE 


5] 


PEASIT BLE “= .<PALSEe 


IF (LON .LE. 217203 FANS SP LON Sere oie THe 


IF (LAT .GE. (33.574 = .0094* TON) EAsteLE ——e Lue 
ELSE DP@LON SCE 25 ie) ee 

IF (LAT .GE. (=7.015 4 .345°DON) SEA SIeee.— — vee 
ELSEIF (LON.2CGE-) 106.2 9s ae 

IF (LAT «GE. (33 ./421)= 20222" Ch} FEAST EGE = Ee 
ELSEIF(LON© <GE™ 106230 ane 

IF (LAT 4GEe 312.742) SPEASTERES= 2 ewer 
BLSELE (iG GE 1 OS .0) a Wet 

IF({LAT .GE. (=31.18 4% .77967DON) SESS Pope woe 
ELSE DPIGON 4GE.) 104 ee ae 

IF (LAT .GE. (=-173.537 54" 129443" 1ON) Sees ahve 
ELSEIF (LON GE) 103 4 aren 

IF (LAT. .GE. (-38.3616 +2. 6500450N) PEAS i Sie ey cela 
ELSEGECLON GE. 103.536) THEN 

LECLAT »GE. 29..0) FEASIBIR = TRUEe 


ELSEIE (LON 2-GE 6102273) aie 


IF (LAT .GE. (166.5411 = 173307*LON) SP EASt eee =e wee 
ELSE DP (GON SCE) 014267 nen 
IF(LAT .GE. 29.837) FEASTSLES= = TRUE. 
ELSE (GON) {GE (2 ee Erie 
IF (LAT .GE. (=/2-6813>4+ 1,009¢7L0ON}) FEAST ELE == sents 
BLUSE IP (DON. GE 6970 464) she 
IF(LAT .GE. (=-104.1222 4 1.3225*LON) ) PEASTBEE = sk Wie 
ELSE TPA GON . GE o920 0 aie 
IF (LAT .GE. (-168.6/709 + 1.9703"L0ON) ) FEASIBLE s— sie wee 
ELSEUE (EON S1GE ws) .lo3 seine 
IF (LAT. .GEe (1.7858 + =324@5*LON)) FEAST BEE = eae 
ENE 
RETURN 
END 
KKKKKKK SMUGSIM SUBROUTINE KEK KKK KEKE KEE KKK KKK KKK KKK KKK KKK RK KKK 
SUBROUTINE SMUGSIM( TEAT, TEON, Terecs ly 
e SMUGGLERS’ SIMULATION - 2/3/87 
: ADAPTED FROM RAND CORPORATION SIMULATION OF VADAP TIVE 
e RESPONSE?) A MODEL OF DRUG INEERDIGRION. 3h eie2. 
C INSTALLATION REMARKS: 
S 
e YOU’LL ALSO NEED TO SELECT "SIZING PARAMETERS” LARGE ENOUCHE 
Cc TO HANDLE YOUR ANALYSIS, BUT NOT SO EARGE AS TO Maree 
CG PROGRAM TOO LARGE TO RUN ON YOUR MACHINE. SEE THE BEGINNING 
CG OF COMMON .FOR™=POR PARAMETERS. 


CERREKEEKKEKRERKEEEKRKEKEKRERERKREEEEKEEKKKEKKEKEKRKKKEKKKEKKKEKEEKKERKKKKKEKEEE 


C PROGRAM VARIABLES: 
INGEGRE COMMON TOR. 


REAL TLAT (MAXASSETS) , TLON (MAXASSETS) 
INTEGER* 40 tk 


CRE KRRRKRKERKRRA RRR R KN AH HR RK DF a ee 


C INITIALIZE THE DATA COLLECTICNPVARRAY Ss 


DO 307 k=]. MAXROUTES 


SV 


oe | 


AMOUNTATTEMPTED(R) 
AMOUNTSUCCEEDED(R} 
PITEMPTS (hae. 
SJCCESSES haa 
DO 20, THISEHASE = 1,MAXPHASES 
PUT EMET Sey Pash kK, THISPHASE) = 0- 
SUCCE Sees PHASE (R,; THISPHASE) = 0. 
20 CONTINUE 
50° CONTINUE 


tot 
oS 


SUCCESSeOsors—0.0 
PrIBORECOSsS=0..0 


SEI SUE Ste TABLE Ore veICHiserOR WPAST SHIPMENTS. 
KMEMORY = ALOG (MEMORYVALUE) / (-LONGPAST) 
BO 40, T=0, LONGPAST 
BRA Bik) Swe AP{ oKMEMORY=REAL (CT) } 
40 CONTINUE 


Coit (hebsovSROUlINE TO COMPUTE ROUTE PROBABILITIES BASED ON 
ASSET LOCATIONS 
CALL PROBCOMP(TLAT, TLON) 


CALL THE SUBROUTINE THA! DOES THE RUN IN. 
CALL PEACE 


DO -GACH YT R TAL: 
DO 50, NTRIAL=1,NUMTRIALS 
Cri. SIMULATE 
ao CONTINUE 


COMPUTE OUTPUT TO ANNEALING ALGORITHM 


Pet iIAb ive rue TOTALS. 
Hora EMe TE D=0.: 

Or Ship PE=0 
LOTATLVEMPTS=0 2 
GOTSUCCESSES=0.. 
mOUSUCCESseost=0.. 

TOIT ALLURECOST=0 . 


SOUEUTE the vLOtAlL NUMBER .OF ATIEMPIS, SUCCESSES, COSTS, (ETC: 
Gres 0 Pekka NUMROULES 
OP ates eas OLA VERVE Ss A TEMP he Gk) 
NOT SUCCESSES = 10 (SUC erssrs + SUCCESSES) 
TOLATTEMPT ED = TOTATIEMPIED. + 
& AMOUNTATTEMPTED (R) 
VOTtSHDP PED =. TOTSHUy PED. + ~AMOUNTSUCCEEDED(R) 
SIUC Ql SSG ISueS SUCCESSCOSTS + 
& PMOUNTSUCCEEDER(R) © “BRUGCOST 
PetiRURECOSTSs= FAILURECOSTS 
( AMOUNTATTEMPTED (R) - 
AMOUNTSUCCEEDED(K)) * 
DRUGCOST 


gad So as Me 


30 CONTINUE 
Pen SUeCESSoCOS P=foT SuCCESSCOs] +SUCCESScCosrs 
FOV AT BURECOST=TOTPAILURECOST+FPAILURECOSTS 


Cela liebe. VERAGE NUMBER OF ATTEMPTS |, Successes, “COsTs, EC. 
Bem TREAL. 
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TOTATTEMPTS .=- TOTATTEMPETS, NUMTR ins 
TOTSUCCESSES = TOTSUCCESSES TUM a e- 
TOTSUCCESSCOST = TOTSUGCESSEGSI, NUiihaike 
TOTFAILURECOST = YOTFAIGUREGOs! NUMER its 
TOTCOSTS = TOTSUCCESSCOsST lc Ve ree 


RETURN 
END 


CRRRERERKEKREKEKEKREREREKKREKRERRKEKEKRERERERREKRKXKKRKKRRAK RK ER ERR KA 5 ee 


CREEK KERKEREEREKRKRERRERKREKREK ERE RR RAR ERKRERR REE RE KER KR Re eee 


SUBROUTINE COMPFACT (RNUMERATOR, RDENOMINATOR, R) 
C GET THE NUMERATOR AND DENOMINATOR OF THE "R" FACTOR, WHICH 
os WILL BE USED TO INCREASE THE PRODAGICD ye eri tp eel lean 
‘S ON ROUTES WITH HIGHER THAN TAVERAGE, TRAV Eare 


CERRRRREKERERRREKERKREREKERERKREKRKEEREEKKERKRRE KEKE EAA ER RRA RO een oe 


C PROGRAM VARIABLES: 


INCLUDE ‘COMMON. FOR’ 
INTEGER* 4 R, RNUMERATOR; ROENOMINATGR Ss 


CRRKKKKKKEKKKKKKEKEKKEKEKKKEKKEKEEKEKKEKEKEKKKEKKEKEKKKKKKKKKKEKKKKEKKEKEKKKKKKKKKKKKKEKKK 


RNUMERATOR = 0 
RDENOMINATOR = 0 
DO LO} S=07 LONGPacT 
RDENOMINATOR = RDENOMINATOR + PASTSHIPMENTS (EK, S) 
10 CONTINUE 
DO 207 -S=0, RECENT Asi 
RNUMERATOR = RNUMERATOR+ PASTSHIPMENTS (RR, s} 
20 CON TINGE 
| RETURN 
END 


CRE RARER EKREREREKRRERRRARAR RRR ARK RK ERR ROR RRR RRR RE eee ee 


CER RAKAKRERREKEKERKEKRKEEKRKEKERKKRKEAEKEKEKEEREEARKEKKEKAKAK ERR A ER eee 


SUBROUTINE GETDATA (GOODINPUT ) 
S READ IN THE INPUT DATA. 


CER E RR ER ERR KER EKER AER RRR EK RR RK RR Re ee 


C PROGRAM VARIABLES: 


INCLUDE ‘COMMON. FOR’ 


REAL Sle Potcey 

INTEGER*4 DAY], DAY2Z7) DAVIN MiIN, PEMe REhe7s 
& P, ROUTENR, 1,0, Pistia eS 
& METHOD, 
& POINT, TY PEAS Se 


CHARACTER*15 POINTNAME (MAXPOINTS) 
LOGICAL GOODINPUD 
LOGICAL TEMPIF 


CEKREKKEREKKEKKREKRKRERKEREREKEREEEKEEKEAEKRKEKKERKKRERKEKKE EERE REKKKKKKKAK 7 


C GOODINPUT WILL INDICATE WHETHER THE DATA WAS CLEAN AND THE 
CG SIMULATION SHOULD BE RUN. 
COCEIMNIPUT == TRUE. 
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OPEN(S5 FibLE=" AVIBAL. BAT”) 


READ( S55, 7} 

READ S59") 

READ ( S>, 4010) STEPSIZEM 

READ Too, 4010) STEPSIZEMIN 

REA 55,4010) STEPS ZELOW 
ev 0’ FORMAT (20X, FS. 1) 

Poon SOr VSrEeP SIZE IN BECREES 
DEtRe— oLEPSrZEM/ 6020 
DEERMING="STlTEPSTZEMIN7 G0: 0 
DELRLOW = STEPSIZELOW/60 .0 


READ (55,4020) BETA 
az 0 FORMAT({20X%, Fs. 2) 


RBAD{(5.,4030)S 
4030 FORMAT (20X, 15) 


READ (55,4040) ALPHA 
4040 FORMAT (20X,F5.3)} 


Ret >5,4050) CUESS 
iO rORMAT (20K,F15.2) 


READ(55,4060) BIGN 
4060 FORMAT (20X,110) 


SEOSE( 55} 


Cc READ IN THE OVERALL SOAR SIMULATION DATA. 
Cah 6, OP ITLE="INPUT DAT’) 
READ(76, *) 
READ (76,5010) NUMMETHODS, NUMPHASES, NUMROUTES, 
& NUMTRIALS, RUNIN, ENDTIME,MEMORYVALUE, IX 
ae FORMAT. (6(15,4X),F10.5, 2X, 15) 
@ CHECK WHETHER TOO MANY METHODS HAVE BEEN REQUESTED. 
TEMPIF = (NUMMETHODS .GT. MAXMETHODS) 
CALL WRTERROR (TEMPIF, GOODINPUT ) 


Cc CHECK WHETHER TOO MANY PHASES HAVE BEEN REQUESTED. 
LEMP E = (NUMPRASES ~Gi > MAXPHASES} 
GALL WRIERROR(TEMPIF, GOODINPUT) 


. CHECK WHETHER TOO MANY ROUTES HAVE BEEN REQUESTED. 
EEMPIPS= ~(NUMROUTES. 2G. MAXROUTES) 
CA WRTERROR(TEMPIT GOOD INPUT) 


Cc CBG Vn HER TOO MANY DAYS HAVE. BEEN REQUESTED. 
TEMP LS = SCE NDTIME. «GT MAXDAYS} 
CALL WRTERROR (TEMPIF,GOODINPUT } 


@ VMeAhVeesURE THE INITIAL SEED IS ODD. (IN CASE WE RUN ON A SUN.} 
ile. — = il), ay aa 
(ale a 


Ss READ IN THE DRUG RELATED DATA RECORDS. 
READ 76, *) 


>») 


READ (76,5020) DRUGNAME, DRUGCCST (EXSHTPMENTINIR Vie 
EXSHI PMENTSIZE, DALEY AMOUNT 


5020 FORMAT “(AlOW4r 07s) 


c 


5030 


READ IN THE METHOD RELATED RECORDS. 
TEMP = NUMMETHODS 
IF (TEMP .GTr. MAXMETHODS) TEMP = MAvMEVHODS 
READ(/6,~) 
DO 74.0. °M= ly Ee 
READ (76,5030) METHODNAME (M), 
RISKCOMP(M) , RISKCOMPEXE (Me CAPCOci hits 
CAPACITY (M) 
FORMATIAS 5X7 4b le) 


40 CONTINUE 


50 


5040 


6016 


60 


70 


* READ 


220 


5080 


IF {NUMMETHODS 2Gi TEME Snel 
DO 50, M=TEMP+1,NUMMETHODS 
READ: (16,5020) 
CONTINUE 
ENDIF 


READ IN THE DAYS WHEN EACH PHASE ENDS? SED UP THES eee 
INDICATING WHICH PHASE IS JN ERPPECT FOKSeA2GCHer A] 
DAY i= 
TEMP — NUMPHASES 
IF (TEMP .GT. MAXPHASES) TEMP — MAX HAsks 
DO ae OF 7 ee 
READS (/ 67,5040) sDAYZ 
FORMAT (15) 
TEMP2Z = SDA v2 
IF (TEMP2 ~.GI. MAXDAYS)> TEMP2 = MAX DAYS 
Te (TEMP2 4 CED Ae a HEN 
WRITE 0/9 25006 aah VEhiE2 


FORMAT (8X, “PHASE * 123)" LASTS VHROUGH DAY rs 
DO 60, DAYN=DAY Ip TEMPEZ 
CURRENTPHASE(DAYN) = N 
CONTINUE 
DAY] = TEMP2 + 1 
ENDIF 
CONTINUE 
CHOSE (6) 
IN THE NETWORK ROUTES AND ROUTE DATA 
OPEN( 96) F1LE= "ROUTES. DAT = 
READ(96G, =) 
READ (96, *) 


DO 250 ROUTENR = 1,NUMROUTES 
DO Ze = i MAX PCinhiie 
NEXT (ROUTENR, §)a— ec 
CONTINUE 
READ (96, 50/75) ROUTENAME (ROUTENR) [ ROUIHCCSt CURE a 
ROUTEMETHOD (ROUTENR) , ROUTELOC (ROUTENR ) 
PORMAT (9A, A107 3X P10. 2 7 2a) 


DOW 2 07 lees 
eM tS) 
COND INGE 
READ(96, S080) FIRST (ROUTENR Presi) —iel 
FORMAT 0-7 bs) 
P=F TRet (ROUTENK) 
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J=1 
240 Ler NES Oo Pen 
NEXT (ROURBERNE FP) = PrS(J) 
P = NEXT CUBOUTENR, P) 
J = eee 
GO-TO - 240 
ENDIF 
250 CONTINUE 


CLOSE(96) 
feReAD IN THE ASSET DATA 


OPEN (2c, PlLE= "ASSET ,DAT”}) 

READ (98,5300)NUMASSETS 

REAU (75,5200 }NUMASSED TY PES 
e200 FORMAT (25x, 15) 


Cc GOUECK WHETHER TOO°MANY ASSETS HAVE BEEN REQUESTED. 
TEMPIF = (NUMASSETS .GT. MAXASSETS) 
CALL WRTERROR (TEMPIF,GOODINPUT) 

C CHECK, WHETRER TOO MANY -AsSEY TYPES. HAVE BEEN REQUESTED. 
MEM -=—(NUMASSETiyYPES .GT. MAAASSEIIY PES) 
CALL WRTERROR (TEMPIF,GOODINPUT) 
BREADS.) 


Dots s0 ASS = Ie aUMASSErS 
READ (98,5320) ASSETNAME (ASSET), 


& LATASSET (ASSET) , LONASSET (ASSET), 

& PSSEVIY PE (ASSET) -ASSETLOC {ASSET } 
220 FORMAT (Go Apa lone (0S) yee iS) 
220 CONTINUE 

READ (9S, *) 

Reet 304") 

BREAD ( 9S, *} 

PEA (oo, } 


DO 3605 TY PE=l) NUMASSEDRTY PES 
READ 93,5325) (RANGE (TYPE, METHOD) , METHOD=1, NUMMETHODS) 
2525 PORMAT (ISA, 318 >. 1,54) ) 
fe) CONTINUE 
READ (9S, *) 
READ( 9S, *} 
READ (IS) ) 
READ (98, *) 
DO 390 TY Pr=1.,.4 
READ (9o,o520) {1 PD (yee, METHOD) MET HOD=1; 3} 
2330 BORMATtTSA, oO (Poss, 24). 
220 CONTINUE 
CEOSE (28) 


eeeeAD IN THE GEOGRAPHICAL DATA 


OPEN 95, LLE="GEOGR.DAT’ ) 

READ (93; *) 
Ber web 95, 5400, END=395) POINT, POINTNAME( POINT}, LAT (POINT), LON( POINT) 
5400 BORMe TO: 1 Sao kA os ov Ee ? soe eS 

COmlOe 372 
pao CONTINUE 

CLOSE (95) 


RETURN 
END 


Sy 


CER EKKERKEKEKEKREKRKEKEREKRKEKERREEKEKEEKRKKEAKKKEKEKEKKEKKARKKKAKKEKEKKKK KE 7 Oe 


CREKKRKKRKKKEKKERKERREKERRKEKRERERKKRERREAR HM RRR OR KK A oe ee A eR eee 


SUBROUTINE SNe iM 
C INITIALIZE FOR THE CURRENT 3SsAyviPhEeroiii 


CEEKKEKKKKKKKKEKKKKKEKKEEKEKEEKKEKKEKKKEKREKKEKKEKKEKKEKEKKKEKKEKKKKKKKKKKKKKKKKKKEKKKK 


C PROGRAM VARIABLES: 


INCLUDE ‘COMMON. FOR’ 
INTEGER = 4. Dik 
REAL RANDNM 


CER KKKKKKKKKEKKKEKKEKKEEKKEKKEEKKEKEKEKEKEKKEEKKKEKKKEKEKKEKKKEKKKKKKKKKKKKKKKKKKKK 


e SET THE*SHIPMENT COUNTER See aro. 
NUMSHIPMENT = 0 


C DETERMINE WHEN THE NEAT SHEPMEND wWikk OGevre 
CALL LRNDPC(IX,RANDNM, 1) 
NEATEVEND = END? IME eo 
NEXTSHIPMENT = - EXSHIPMENTINTRVL * ALOG(RANDNM) 
IF (NEXTSHIPMENT .b7. NEXTEVENMI seinen 
NEXTEVENE = NESE So eve bn 
ENDIF 


C GET THE PASTY FROM THE LONG GRUN UN: 
DO 30, BD=0;, LONGPAST 
DO 20, R=1,NUMROUTES 
PASTSHIPMENTS (hb) = PEACE SH LEVEN TS) 
PASTFPAILURES( R72) e=s PEACE PAIMURES (hb) 


208 CONT INUE 
30 CONTINUE 
RETURN 
END 


CRERRKKKRKKKEKKRK HERRERA KERR RRR RR RE RK RRR AR RO RRR ee eee 


CERRRKREKKKEKKKEKKEKKEKKEEKEEKKEERERKKKEKKKEEKKEEKKEEKKEKREEKREKKREAKA AX A A 


SUBROUTINE PEACE 
C DO THE INTIAL RUN@ Ls 


CERKKRKKKKEKKKKKEKKKEKKKKKEERKEKKKKKKKKKKKAEKKEKKKKKARKEKKKKKKKKKKKAK KX 7 ee 


C PROGRAM VARIABLES: 


[NELUDE. COMMON: FOR} 
REAL AMOUNT, PROBCAUGHT, RFACTOR, TEMPTIME, RANDNM, ACTUALRISKCOMP 
INTEGER*4 D,R, RDENOMINATOR, RNUMERATOR, T, TRIPS 
LOGICAL ER DEeSUCeESs 


CREKKKKKKKEKKKKKKKKEKEKKKEEKKKKKKKKEKEKEKKEKKEKEKKKEKEKKKKEKKKKKKKKKKKKEKEKKKKKKKKKKEK 


C DETERMINE, THE TIMES Or Sink shi Ses per Eee 
NEAXTEVENT =— RUNING 130 
CALL LRNDPC (IX, RANDNM, 1) 
NEATSHIPMENT = -EXSHIPMENTINGRY ibe VALOG(RANDNM) 
IF (NEATSHIUPMENT iN VED) Ebish| 
NES TEVEN S=— NEXT Sree rinnit 
ENDIEE 
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10 


20 


aH 


PNITTALIZE THE “ARKKAYS DESCRIBING THE=PAST . 
DO 230, D=0, LONGPAST 
DO 2207 ke OMROUTES 
PASTCHUPRMEBENTS{R,D) = 0 
PAS WnUR ES (kh, Diy = 0 
CONTINGE 
CONTINUE 


FORK -EAGH DAY OP] RUNIN =. 
DO 100, DAYNOW=1,RUNIN 
TiS Pras ta 1 


Silt THES ARRAY S DESCRIBING. THE PAST. 
DO=Z207 R= NUMROUT ES 
DOW LU D=LONGPAS iy - 1 


Pas sive Mblile (hp) = PASTSHMPMENTS(k,D-—1} 
Pesta bURESth, Oo) =) PAStrPALLURES(R,b=-1} 
CONTINUE 
PASTS EEMENTS (0) = —0 
Peo eens (ik .0) 52.0 


CONTINUE 


EOOP THROUGH THE DAYS SHIPMENTS. 
TPN UE eT. 2 GE a MoAYNOW 4-120) 9GO° TO, 100 


Cee TiEr OU Y OF SUE SNE xd eset EMenN Ty. 
AMOUNT = EXSHIPMENTSIZE 


SEPEEC Telit thOU TEs OseBe USED. 
CALL SEDROUTE (RK, PRIPS; AMOUNT ACTUALKISKCOMP) 


COMPUTE THE NUMERATOR AND DENOMINATOR OF THE 
Uh T EAR 
CALL COMPFACT (RNUMERATOR, RDENOMINATOR, R) 


HORS Mreh ether Ik BOULRED TO GE AMOUNT SHEP PED... 
DO 07 el =e eS 


COMPO wis Kh =PACTOR:. 
RPACL ORs ——20 
IP“ CREENOMINATOR. -<C)..0.0). THEN 
RFACTOR = RNUMERATOR/RDENOMINATOR 
TEE TIME = Eee vei 
Ie ATEHPIIME 26). RECENT PASI)? “Hen 
PS VTEMP IME =.GT 28 LONGPAST) 
TEMPL EME." LONGCPAST 
REACTOR —~ erAClTOR =~ “TEMPTIME/RECENTPAST 
ENDL 
ENDIF 
iPr PACTOR. it 0) RF ACIORS=, 1.20 


COMPUTE =e PROBABILITY OF INTERDICTION. 
IP OCP ROEMVIER@ICT (KR) THiPSPHASs )-Gr..9939). THEN 
PROBCAUGHT = 1.0 
SESS. 
PROBCAUGHT (=. 15 05— EPXP( REACTOR -ALOGC1...0: = 
PROBINITERDICT(R, THISPHASE})) 
ENDIF 


DETERMINE WHETHER THE. TRIP wWASsSUCCESSFUL. 
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CALL LRNDPC (IX, RANDNM, 1) 


TRIPSUCCESS = (RANDNM GE Sy ROpGAUCn 
GC DO THE BOOKKEEPING. (NOT VER Yer ENS? VE, 
C DURING THE RUNS. 
PASTSHIPMENTS (R70) 9=) Asi SH TEMES ea eee 
IF (.NOT UR resveere 
& PASTFAILURES(R,0)) =, PASTFAILURES (RAG oe 


RNUMERATOR = RNUMERATOR + 1 
RDENOMINATOR = RDENOMINATOR + 1 
70 CONTINUE 


C GET THE TIME AND Ty PE SOT rb Ue Xai ra lest teshaplee 
CALL LRNDPC( 1X, RANDNM, 7) 
NEATSHIPMENT = NESTE ee 
& EXSHIPMENTINTRVL * ALOG (RANDNM) 
NEXTEVENT = RUNIN + 1.0 
IF (NEXTSHIPMENT “(int NEXTEVENT ewe 
NEXTEVENT = NEXTSHIPMENT 
ENDIF 
GO; TO 20 
100 CONTINUE 


S SAVE THE LAST LONGPAST DAYS FOR USE INITIALIZING EACH TRIAGE 
DO 130, D=0, LONGPAST 
DO 1207) Rai NUMROUTES 
PEACESHIPMENTS (R, BD) =PASTSHIPMENTS (Rep) 
PEACEPAILURES (RD) PAST PAT RURES (hye) 


120 CONTINUE 
130 CONT INUE 
RETURN 
END 


CERRKKKKKKKEKEEKEEEEREKKREEEREKEKEERKRKEKREKAEKEKREAKERRERAEKEKEKKKEEX EE ASE 


CERRRKKKKKRKREREKKEKEKRKRERREAKEKRERREREKRERAEKREKREAKKEKRKKEKKAKERKKAAKAEKEKRRAKKEAX A 7 ee 


SUBROUTINE SELROUTE (RCHOSEN, TRIPS, AMOUNT, ACTUALRISKCOMP) 
C SELECT THE ROUTE PORK THE NEAT snr eienie 


CREE KKK EKER KKK EEK KEKE KKK KKK KKK KEKE KE KEKKKKEKAKKKKKKEKEAKKEKKKEKKKEAKKKEKKKKEKAEKKEKKKEKEK 


C PROGRAM VARIABLES: 


INCLUDE, ~COMMON FOR: 
REAL AMOUNT, CUMPROB, ROUTEPROB, COSTROUTE (MAXROUTES ) , 
ACTUALRISKCOMP, PROBCAUGHT (MAXROUTES) ; 
TEMPNUMERATOR (MAXROUTES) , TEMPDENOMINATOR (MAXROUTES ) , 
& TOTCOST , WEIGHT, TEMPRISKCOMP (MAXROUTES) 
INTEGER ~4 “KCHOSEN, TRIPS, 5METIOD n> 


6 Sh S| 


CEE KKKKKEKEKKEKKKEKEKEREKKKKKEKEKEKKEKEKKEEEKEKKEKEKEKKEKEKEKKKEKEEEKKEKAKKEKEKAEKEKKEKEKES 


c INITIALIZE THE NUMERATOR AND DENOMINATOR FOR EACH ROUTE. 
DO] 07 k=), NUMECUTES 
TEMPNUMERATOR (5) =~ 0. 
TEMPDENOMINATOR(R) = 
COSTROUTE (hw —- 70,0 
10 CONTINUE 


& FOR EACH DAY TO BE CONSIDERED, SACD diets CONTRIBUTION: 
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SES Reue na 


(ee) 


20 


& 


& 


BOS. .S=0 /LONGPAST 


WEIGHT = BAPTABLE(S) 
DO 20, R=1,NUMROUTES 
TEMPDENOMINATOR(R) = 
TPEMERENOMENATOR (8) -+ WEIGHT’ PASTSHIPMENTS (k,S} 
TEMPNUMERATOR({(R) = TEMPNUMERATOR(R) + 
WEIGHI™ PAST AaGURES (Re S) 
CONTINUE 


eS CONTINUE 
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PORs PACH SOULE 2 COMPUTER THE PERCEIVED PROBAB Iba TY OF sBE TING 
CAPTURED ANDORENCEH ATE (bP XPPCLE Ds COst sor USING THE ROUTE. 
THES PROBABIIGITY A ROUTE WILL BE CHOSE WikLs BE PROPORTIONAL 
TO STH SVE Ret bree EDC OS OF USiniG THE ROUTE. SO. GET THE 
TOTAE Ory THE siNVER SE AEX PECTED COSTS OF USTNG: BACH ROUTE . 
POTGOST = “U0 
DO R=) .NUMROUTES 

RMBTHOD = ROUTEME THOR (R) 

TRIPS = iNT AMOUNT CAPACITY (RMETHOD) ~+..0..999 ) 

IF (TEMPDENOMINATOR(R) GT. 0.001) THEN 


PROBCAUGHT(R) = TEMPNUMERATOR(R) /TEMPDENOMINATOR (R) 
Elbor, 

PROBCAUGHI(R}) =) 00 
ENDL 


iP et RUSK COMPEAP{EMETHOD): 2UE. 020001) “THEN 
TEMPRISKCOMP(R) =RISKCOMP (RMETHOD) 
EUs: 
TEMPRISKCOMP(R)=RISKCOMP(RMETHOD) * 
Gee PROBG-UGH IT GR RISRCOMPEXAP(RMETHOD)) 
ENDIF 
COSTROUTE.UR) .= 12077 URS] (PROBCAUGCHIAGR = CAPGOLT( RMETHOD) + 
ROUTECOST (kK) es TEMPRISKCOMP(R) } + 
PROBCAUGHT (R) * AMOUNT * DRUGCOST ) 
TOTEGOST = Telecel 2 ces TRourte rR) 
CONTINUE 


NOW CHOOSE THE ROUTE, WHERE THE PROBABILITY OF CHOOSING: A 
GIVEN ROUTE -15~PROPOKTIONAL TO THE. INVERSE COST OF USING 
THAT ROUTE. 

CALGACRNBES( 14, ROUTEEROB 31} 


CUMPROB =-0.0 


= 


EO CONTINUE 


eae. ae | 
CUMPROB = CUMPROB + COSTROUTE(R) /TOTCOST 
(CGUMPROB. JLT. ROUTEPROB) GO TO 50 


RCHOSEN = R 


RMEBUGHOD = — ROUTEMETHOD(RCHOSEN} 
TRIUPS, =sTNt (AMOUNT, CAPACITY (RMETHOD) + 0.999) 
AC TUALRITSKCOMP=1TEMPR LSKCOMP( REHOSEN) 


RETURN 
END 


eee RK ARR KARR RK RRR RA eA ee REE EE REA RRR ERX 


Sate KKK KKK KARE RKAKKKE KA AKEAKEKRE KARR EKKRRKAKRREKKKERKEKKEKKEKEKKKEREREE 


aA 


SUBROUTINE SIMULATE 
THIS SUBROUTINE PERFORMS THE SIMULATION FOR A SINGLE SAMPLE 


POINT: 
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CEE REKRKKEKEKEKRKEKKKKEKKKEKREKREKKKEKKKKKREKRAEEKRAKKKKEKER RK KEK KKK KKK KK 5 Ree 


C PROGRAM VARIABLES: 


INCLUDE .* COMMON ROR: 
REAL AMOUNT, PROBCAUGHT, RFACTOR, TEMPTIME, RANDNM, ACTUALRISKCOMP, 


& TIMESHIPPED(MAXSH IT PMENTS) 
INTEGER*4 D,R, RDENOMINATOR, RiUMERATOR Tees, 
& ROUTEUSED (MAXSHIPMENTS) 


LOGICAL TRIPSUCCESS, SUCCESS (Mast ME iis: 


CEE RKEKK KKK KEKE KEK KEK KEKKEKKKEKEKEKEKEKKEKEKEKEKEKKKKKKKKKKKKEKKKKKKK KK KK KKKKKKAKKKKK 


@ 


CeO) ct eG) 


INITIALIZES OR Thise Rial 
CALL INITS iM 


FOR BACH PAY. 
DO 100) DAYNOW=-0) ENDEIME 


GET THE POINTER INTO THE ARRAY OF INTERDICTION PROBABIL Agree 
THISPHASE = CURRENT PHASE (DAYNOW) 


SHIFT THE PAST HistTonyY ARRArsS. 
DO 720, (R=) NUReurE. 
DO 10>, D=LONGPAST) 1-1 


PASTSHIPMENTS (RK, D) = PASTSHITPMENTS (hyp!) 
PASTFAILURES(R,D) = PASTFAILURES(R,D-1) 
10 CONT INUE 
PASTSHTEMENTS (he U aan! 
RBASTFAITLURES( RoC) ae 
20 CONTINUE 


LOOP THROUGH THE EVENTS tiad. HAPPENS TORAY = 
oi IF (NEXTEVENT .GE. DAYNOW + 1.0) GO TO 100 


DETERMINE THE NEAT TYPE OF DRUG TOC BE VSI PrED eA 
HOW MUCH. 
AMOUNT = EXSHIPMENTSIZE 


SELECH it OUTE. 
CALL SELROUTE(R, TRIPS, AMOUNT, ACT UAER TS) GoM) 


GET THE NUMERATOR AND DENOMINATOR OF THE “R" FACTOR. 
CALL COMPFACT (RNUMERATOR, RDENOMINATOR, R) 


MULTIPLE TRIPS (MAY BE REQUIRED TO Stra tieeoly Ey 
AMOUNT OF DRUG ON THE SELECTEDSROUTE ) GEO 
CAPACITY CONSTRAINTS. THE VARUAEEE TRiPss2 iwi 
THE NUMBER OF TRIPS THAT WILDL PEE REOU [RED bo 
THAT THIS VARIABLE WILL USUALLY BE CONTAIN A ONE 
AND THAT AS A RESULT DO-LOOP 70 WILL USUALLY BE 
EAECUTED ONE ONCE. 

DO. 70. T= hes 


COMPUTE ThE sek. AChOke 
REACTOR. — 1.0 
TF <(RDENOMINATOR 3.Gl) 020 eV neEn 
RFACTOR = RNUMERATOR/RDENOMINATOR 
TEMP IME S=. NE eke a 
IF ({(TEMPTIME, GT >] RECENTPAST) Sha 
TF ClrENPT IME Gi LONCrA 
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70 


100 


GO T 
COMNLINUE 


TEMPTIME = LONGPAST 
REACTORS = REACTOR * TEMPLIME/RECENT PAST 
END EE 
ENDIE 
ieee lO bid OO eRRACTOR .=*.1..°0 


COMPOMESITHE PROBABIL TY .Or- INTERDICTION. 
COMPUTES iE Se hOBAG Liiiiy sor? INTERDICTION. 
he POR eRe or eR PH SPHASE) 4GE779999)— THEN 
PROEPCAUGHI »=—120 
EGSE 
PROBCAUGHIA = 10 Vex (RRACTOR? AlLOGt 1. 0 = 
PROBINTERDICT (8, 1HISPHASE )})) 
ENGEL 


DETERMINE SWE fate REP WAS SUCCESSFUL), 
CALL LRNDPC (IX, RANDNM, 1) 
TRIPSUCCESS =) (RANDIM ~ GE. )PROBCAUGHT } 


DO THE wREQUIRED VEOORVEEPING. 
NUMSHIPMENT = NUMSHIPMENT + 1 
SUCCESS. NUMenI PMENT) = TREPSUCCESS 
TIMESHTEPPED (NUMSHLEMENT ) 2= NEATEVENT 
ROUTEUSED( NUMSHTPMENT a == Rk 
PEEVE TS (Kk) =a treMePTsS (Kise 1.0 
AMOUNTATTEMPTED(R) = 
AMOUNTATTEMPTED(R) + AMOUNT/TRIPS 
ATTEMPT SBY PHASE (Rtn SPRASE) = 
Po. TEMPTS BY PHASE Gh, VHUSPHASE) +1 
PEA TRIP SUCeCESe je THEN 
SUCCESoCCols=SUCCESSCCsST SS ACTUALRITSKCOMPT 
ROUTEGOSTAR) 
SUCCESSES (Rh) 2= SUCCESSES (RK) y4 120 
AMOUNTSUCCEEDED(R) = 
AMOUNT SUCGERDED(R) “+ “AMOUNT /TRIPS 
SUCCESSES SB VPs oe tn, IHiESPHASE A= 
SUCCESSES RY PHASE (Re THISEHASE): +. 1 
ELSE 
PAStPAtEURr Ss hoi eP ASA LUE ES tie 0) 4 ol 
FATLURECOSTS=FAI LURECOSTS+ACTUALRISKCOMP+ 
ROUTECOST UR) 4 
CAPCOST -ROUTEMETHOD (R})) 
ENDIF 
PAST SHITPMENTS (R70). ="“PASTSHEPMENTS (RR, 0) + 1 
RNUMERATOR = RNUMERATOR + 1 
RDENOMINATOR = RDENOMINATOR + 1 
CONT INUE 


DETERMINE SWHEN athe. NEAT SHEPMENT Or THIS DRUG WILL BE. 
CALL LRNDPC(IX,RANDNM, 1) 
NEAPSHLeMENT = NEAT EVENT = 

EXSHIPMENTINTRVL * ALOG (RANDNM) 


DETERMINE WHAT THE NEAT SHIPMENT WiLL BE. 

NEATEVENT -=cENDT IME 422130 

DF (NEAT SH PeMENT= Lis NEATEVENT): “THEN 
NEATEVENT -="NEATSHEPMENT 

ENDIF 

© 3.0 
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RETORN 
END 


CE RRR ERR RR BK RR KR eee 


CR RRR RRR KK KER REE RRR KR eee 


SUBROUTINE WRTERROR (ERRORCONDITION, GOODINPUT) 
‘S THIS SUBROUTINE WRITES Se ee ere MbSeeCe ir NEEPER- 


* 


CEKKKKKEKEKKEKEKREEKERKKKEKEREKREKERERKERRAKKKKAKKEKEKEKRK KKK KKK KX KK KOK XX eee 


C PROGRAM VARIABLES: 
LOGICAL ERRORCONDITION, GOODINPUT 


CERKRKEKEKKEREKRKEKEKKEEEKEKREEKKKKERKKREEKKKREKKEKKKEKKKKKEKKEKKEKAK KE 5 EX Oe 


IF (ERRORCONDITION) THEN 
WRITE? 76000) 


6001 FORMAT “(9 *4BRROR sce) 
GOCDINEUT =s2.25sSre 
BIND ITE 
RETURN 
END 


CEE KKKKKKRKEKKEKKEKKKEKKEKEKKKKKKKKKKEKEKEKEKKKEKEEKKEKEKKKKKEKKKEKEKEKKEKKKKKEKEKKEKEKE 


KEK KKK KEKE KKK KEKE KEKE KEKE KK KEKE KK KEKE KKK KKK KKK KEKE KKK KKK KEKE KEKE KKK KKK KEKE KKK KKK KEKE EK 


SUBROUTINE PROBCOMP(TLAT, 1hen) 
THIS SUBROUTINE, COMPUTES PROBAR TA iissore lie Olt Tone 
ROUTES AS A FUNCTION OF ASSET spiel Nero ROUlE: 
PROBABILITY OF INTERDIGT ION Ts) teu te stiles ii. 


+ + + + 


INCLUDE “COMMON TiCke 


REAL M1,M2,3B1, B2, TEMPP], EATINT SUONINT DST] DEND, BENGI AH Piste 
& TLAT (MAXASSEDS); 1 LONI ASsokl Ss) Coble, FCcrr 
INTEGER SE) Orr 


* COMPUTE PROBABIEI IY OF INTEREITCT TON PORE ereenhOU le 


DO 405° R=1) NUMRCUTES 
DO 400 NSE NUMPRASES 
PROBINTERDIG) (a) eee 
400 CONTINUE 
405 CONTINUE 


* FOR EACH ROUTE COMPUTE EQUATIONS FOR ROUTE Secvria= 
DO 440 R=1,NUMROUTES 


P=FIRST(R) 
O=NEXT (RR, FP) 
410 CONTINUE 
IF ( LONCP) SsNE EON:(o) . AND. 
& LAT (P) .NED GAY? (O)) THEN 
Mi= (LAT (Pelee EG) aye 
& (LON(P) - LON(Q)) 
M2=-1.0/M1 
BI=LAT (CP) = ie bone? 
ENDIF 
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TEMPER =O 0 
DO 420 ASSET=1, NUMASSETS 


TE EON CE) SNE. LON (0) SSIES 
& Bee?) NE bel (0) THEN 
Bea SSE) (ASSET) - v2 TONASSE!] (ASSET) 


Eat Tie M2 82 Mie) (M2 -M 1) 
IO) SO = OsME ee gu ull 
EDSE Tit hel Cr) SEO LON GC )2 (THEN 
LATINT=LATASSET (ASSET) 
LONINT=LON (P) 
ESE 
DATUNT=LAT (Pt 
LONINT=LONASSET (ASSET) 
ENDL 


fe OETERMINE CLOSEST POINT BETWEEN ASSET AND ROUTE SEGMENT 
Det —=GSbis Loni r)] barr), 


& LONINT, LATINT) 
DEND=GSDIS (LON(Q) , LAT (Q), 

& LONINT, LATINT ) 
EENGTINeCSPiS (EON TP). bare), 

& LON(O), LAT (Q) ) 


TP VbSt SLE. LENG THe ab WEND 7. bE? aLENG TD THEN 
DIST=GSD1I5{ LONASSET (ASSET) GATASSET( ASSET), 
& LONINT, LATINT ) 
BESE LE Qest «iit .) END). THe 
Dist—Cobis (LONASSED (ASSET); LATASSET (ASSET), 


& DOM) AT Ce) 
BESE 
Dist-GsSbisc | LONASSET( ASSED) (LATASSET (ASSET), 
& LON (Q) , LAT (Q) ) 
ENDIF 


fee corUlE PROBABILITIES (LNDEPENDENCE ASSUMEDO 
iPOD iS apo 
& RANGE (ASSETTY PE( ASSET) , ROUTEMETHOD(R))) THEN 
A= 0 ee (ASSEVIT YEE CS ce OUT EMETHOLYR)) 


ELS& 
TP (RANCH (ASSET IEE (issb? es ROUTEMS Ino Uh) )- 2bOr 0.0) 
& RANGE (ASSET YPE(ASSED) ROUTLEMETHOD(E)} = 0.0001 
PCTR= (RANGE (AsSSEtIYER ASS Rou TEM TOD Ri g=DIsST) / 
& RANGE (ASSEITY PE ASsEr) ROUTEMETHOD( R}) 
AJ = 1205 =— TPD( ASSENT VY PE VAS oe!) ROUTEMETHOD(R) ) 
& nap CE Ci) 
ENDIF 
AL yale Oe eee et Op ec EMEP Iya) 
420 CONTINUE 
DO’ 430 N=), NUMPHASES 
PROB UNTER DIGI CR AN) 2] 100 = ole = TEMPP 1) 
& (22 0°= PROBUNTERDIET(R,W).) 
430 CONTINUE 
J 1G, 


Ot= 7 hE es 
te NE. 0) 9 GO sTO741¢G 
440 CONTINUE 
RETURN 
END 


Saree RAK RAK KAR XKKAARAEAKRRRAKRERERAREREKAREE ER EKER KAKKAKEREREKREREREEER 


Peer ONeGIMION CGSDIS{ xl, Yl, 42, Y2 } 
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* - FUNCTION COMPUTES THE GREAT -CIRCLE DISTANCE ieee 
* BETWEEN TWO POINTS EXPRESSED IN» GAT ITUDE SAND eioOner Tyee 


REAL*4 Kale eee 
REAL*4 RARC, RLAT1, RLAT2, RLON1, RLONz, RPD 


*-—--LOCAL DATA INITIAbI ZA Tron 
DATA RPD /0.1745379433017 307 ea0e 


*=—-GREAT -CARCLE SOUPS TANCE 
* 
m INPUTS ARE TWO POINTS EXPRESSED AS (Wo LONGITUDE, N22 bATIIGpas 
*=-= CONVERT DECKREES, TO KAUIANS 
REON ==. hPa 
REA =e ea 
RLON 2 = RP De we 
RLATZ = Reps wZ 
*--—-SPHERICAL ARC LENGTHS IN RABIANS 
ARGUMENT = (SIN(RLAT1) *SIN(RLAT2) 
& +COS (RLAT 1) *COS (RAT Z) “COST REONT=RLONZ), 


IF (ARGUMENT. .GT..1.0) ARGUMENT =a ec 
LFP( ARGUMENT 27. —1 0) SAR CUME Nao 
RARC = ACOS{ARCUMENT) 


*---ARC LENGTH IN STATUTE MILES WITH ROpeiReCORRECTION FOR VOELATE I= 
. OBLATENESS CONTRIBUTESSEESS (HAhees 25s CoRREC IGN: 

i IN U.S., DARC IS APPROAIMATI ieee URAC 
* 


GSDIS* = 395770807 Rane 


RETURN 
END 
oy Poy ee RANDOM NUMBER GENERATING SUBROUTINES “7 ******s****** 7) 


ADAPTED FROM LEWIS, ORAV,AND URIBE ENHANCED SIMULATION AND 


STATISTICS “PACKAGES [Ret 27] 
THIS PROGRAM WILL GENERATE A VECTOR OF NORMAL RANDOM VARIABLES 


ACCORDING TO- THE SsiINE-CesS tik ME Tne. 


CuSrerene 


SUBROUTINE GHOREC (TSeEbpy4 NN} 
INTEGER. WN iin Seer 

REAL) ACN} 7 U(2) pSV ST Ak 
DOUBLE -EREGISION 7 

DATA Bly 2 14 2658-67 ey 


IND Sse 
DO; 100 ei ia 
JUNIE ee NIG 
[TF UND eCE SGCOlre.ce 
10 CALE VLENDEC (TShED 072, O70) 


S =S0ORTM-2 7 ALOG (Uta 
WS 2 Pl) 
XSTAR = S*COS(W) 


Atl) = S*sSinwy) 
GOTO 100 
20 A(I) = ASTAR 
100 CONT LNW 
RETURN 
END 


66 


SUBROUTINE —LRNDES wel seeb UN) 


INTEGER Nee Pore 

REAL U(N) 

DOUBLE PRECISION boii. DSEED, <De1 
CG Bead 2 Sl 
‘e Pete =e oo 


BATA DSIMiL7 2A a 33 64720, 
WATA 0S 1) 72a 4es Gab, 
PSEED == UE EE (1 CEED) 


Demo et i 
& Pst = Mee s507 Ves? 2b) Shieh ms iMi) 
DSEED = DMOn (les Wj. Due bSEED, D2iM1) 
UC sD SEED7 bo 
Ss CONTINUE 
PEED = INE (PSEED> 
RETURN 
END 
tala COMMON . F'OR 
E Sod “PAMETERS: 
INTEGER* 4 LONGPAST., RECENTPAST ,. MAADAYS, 
& MAXMETHODS, MAXPHASES, MAXROUTES, MAXSHIPMENTS, 
& MAAASSETS, MAXASSETTY PES, MAXPOINTS, 
MAXPTS 
PARAMETER (LONGPAST=120 ,SRBeCENTPAST=20, 


MAXDAYS=730, 

MAXMETHODS=3, MAXPHASES=12, MAXROUTES=100, 
MAXSHIPMENTS=6000,MAXASSETS=100, 
MAXASSEITY PES=4, MAXPOINTS=50, MAXPTS=8) 


gl a a el 


S COMMON VARIABLES: 


c REAL VARIABLES: 

REAL AMOUNTATTEMPTED (MAXROUTES) , 
AMOUNTSUCCEEDED (MAXROUTES), 
ATTEMPTS (MAAROUTES) , 
ATVTEMPISE VPHASE( MAX ROULES  VAAPHASES) , 
CAPCOST (MAXMETHODS ) , 
CAPACITY (MAXMETHODS) , 
DAILYAMOUNT, 
DRUGCOST, EXPTABLE(0:LONGPAST) , 
EXSHIPMENTINTRVL, 
EXSHEPMEN Rad Zi, 
FAILURECOSTS, KMEMORY , MEMORYVALUE, 
NEATEVEND NEATSHIeEMENT, 
PROBINTERDICT (MAXROUTES,MAXMETHODS) , 
RISKCOMP (MAXMETHODS) , 
RISKCOMPEXP (MAXMETHODS) , 
ROUTECOST (MAXROUTES) , 
SUCCESSCOSIS, 
SUCCESSES (MAXROUTES};; 


2G eG hg a a Soa a Soe a a a 2 eh a ag fee a Od 
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Seg Ie Sto JN oe Melk Sn Fe So 


SUCCESSESBY PHASE (MAD ROUTES, Mae esr s)] 

LAT (MAXPOINTS) , LON(MAXPOINTS), 

LATASSET (MAXASSETS) , LONASSET (MAXASSETS), 
DELR, DELRMIN, DELRLOW, 

BETA, ALPHA, RANGE (MAXASSETTY PES , MAXMETHODS) , 
IPD(MAXASSEDTYPES Mein HOS) SGubSs 


INTEGER VARIABLES: 


INTEGER*4 


m 


Br Re (Ro OR? Bre ye ea. Bag 


CURRENT PHASE (0 :MAXDAYS) , DAYNOW, ENDTIME, NTRIAL, 
NEXTEVENTTY PE, NUMMETHODS , NUMPHASES, 
NUMROUTES , NUMSHIPMENT, NUMTRIALS, 
PASTSHIPMENTS (MAXROUTES, 0: LONGPAST); 
PASTEAILURES (MAXROUTES 0 -LeNerast © 
ROUTEMETHOD (MAXROUTES) , RUNIN, 

THIS PHASE: 
PEACESHIPMENTS (MAXROUT ES, 0 - LONGEASTy 
PEACEFAILURES (MAXROUTES , 0: LONGPAST )7 
NUMASSETS, ROUTELOC (MAXROUTES), 

ASSET LOC (MAR ASSETS) =, 
ASSETTYPE(MAXASSEITS) SNUMASS EA Ey PEs, 
S, NEXT (MAXROUTES ,MAXPOINTS) , 

FIRST (MAXERCUIES } 7 BIGN 


INTEGER VARIABLES: 


TN ECER 4 


Ie Sieve] 6, 


CHARACTER VARIABLES: 


CHARACILER © 4.0 


CHARACTER 212 
CHARACPER SiS 


METHODNAME (MAXMETHODS ) 


ROUTENAME (MAXROUTES ) 
ASSETNAME (MAXASSETS ) 


COMMON AMOUNTATTEMBBED, AMOUNTSUCCEEDED ATTEMPTS; 


Sl Sa N tal aa gM) Sg oe NE oe Nal og NO op og Wi ae a ag VE og al ee eg a Loe td Mac a Sag a Se I a sg Sig Seg og. 


ATTEMPTSEY PHASE, 

CAPCOSE, 

CAPACITY, 

DAILYAMOUNT, 

PRUGCOST VEAP TABIE, 
EXSHIPMENTINTRVL, 

EZ SHLEMENTS 128; 
FAILURECOSTS, KMEMORY , MEMORYVALUE, 
NEATEVENT, NEATSHIPMENT, 
PROSIMTERE £e ie 

RISKCOME, 

RISKCOMPEXE, 

ROUTECOS 

SUCCESSeOs ls) 

SUCCESSES. 

SUCCESSESBY Pie ee, 

LAT, LON, 

LATASSET, LONASSET, 

DELR, DELRMIN, DELRLOW, 

BETA, ALPHA; RANGE, LED, CUBS. 
CURRENT PHASE, DAYNOW, ENDT IME, NTRIAL, 
NEATEVENTIY PE, NUMIETHODS NUMPHASES, 
NUMROUTES, NUMSHIPMENT, NUMTRIALS, 
PASTSHIEMENTS. 

PASTFAILURES, 

ROUTEMETHOD, RUNIN, 

THISPHASE., 
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Ra SS ee ey a 


BEACESHIPMENTS, 

PEACEE AILURES, 
NUMASSETS, ROUTELOC, 
ASSETLOC, 

ASS iy PES NOUMASSH iy PES, 
Sibeer, | LKST, BIG, 
oD, 

METHODNAME, 

ROUTENAME, ASSETNAME 
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